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Abstract 

Applications imposing SU(3) chiral symmetry on non-relativistic field 
theories are considered. The first example is a calculation of the self- 
energy shifts of the spin-| decuplet baryons in nuclear matter, from the 
chiral effective Lagrangian coupling octet and decuplet baryon fields. 
Special attention is paid to the self-energy of the A baryon near the 
saturation density of nuclear matter. We find contributions to the mass 
shifts from contact terms in the effective Lagrangian with coefficients 
of unknown value. As a second application, we formulate an effective 
field theory with manifest SU(2) chiral symmetry for the interactions 
of K and r\ mesons with pions at low energy. SU(3) chiral symmetry 
is imposed on the effective field theory by a matching calculation onto 
three-flavor chiral perturbation theory. The effective Lagrangian for the 
7rK and tti] sectors is worked out to order Q 4 ; the effective Lagrangian 
for the KK sector is worked out to order Q 2 with contact interactions 
to order Q A . As an application of the method, we calculate the KK 
s-wave scattering phase shift at leading order and compare with chiral 
perturbation theory. We conclude with a discussion of the limitations of 
the approach and propose new directions for work where the matching 
calculation may be useful. 
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Chapter 1 
Introduction 

I once heard a wise man say "In the history of scientific endeavor, no problem has con- 
sumed as much of mankind's resources as the understanding of nuclear forces," or some- 
thing like that. In any case, for all the considerable effort poured into solving the mysteries 
of the strong interaction, a number of significant problems remain. The fundamental the- 
ory of the strong interactions, Quantum Chromodynamics (QCD), is solidly established 
as a pillar of the Standard Model of particle physics and, to the extent that QCD is a 
renormalizable gauge theory, is well understood in the perturbative regime. As a theory 
with asymptotic freedom [I], ^J, QCD is perturbative in the high-energy regime; for low 
energies the coupling constant of the theory becomes large and perturbative treatments 
break down. Some features of non-perturbative QCD which are still not fully understood 
are quark structure of hadrons, dynamical symmetry breaking, and quark confinement. 

The special difficulties accompanying the non-perturbative regime require special meth- 
ods for working within the theory. One direct approach is to formulate QCD on a lattice 
of space-time points and use numerical techniques to perform the functional integrals. 
The lattice QCD method || || has great potential for shedding light on many of the 
unanswered questions of QCD and has become an industry unto itself. A complementary 
approach is to focus on the long-distance physics, and base the field theory description of 
the physics on the directly observed degrees of freedom. Approaches of the second type 
are generically called effective field theories and rely on a two-part foundation. I refer to 
the first important concept as the Weinberg Hypothesis ||: the only content of quantum 
field theory (apart from the choice of degrees of freedom) is analyticity, unitarity, cluster 
decomposition, and the assumed symmetry principles. As a consequence, we can describe 
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the strong interactions in terms of hadron degrees of freedom provided we write down 
the most general Lagrangian consistent with the symmetries of QCD. The second key 
concept is that we must identify some expansion parameter, typically a small momentum 
or energy scale, which permits us to calculate to any given order in the expansion with 
a finite amount of work. The predictive power of an effective field theory arises from the 
combination of the two underlying concepts; the symmetry restricts the parameters of the 
theory to a meaningful set and the expansion parameter allows a systematic framework 
in which we can include all contributions of a given order and estimate the size of the 
contributions we have neglected. For a more detailed discussion of effective field theory 
in general, see references || |7|, |8| . 

In the case of QCD, the hope is that by experimental determination of the param- 
eters of the effective Lagrangian, we can learn about the underlying theory in the non- 
perturbative regime, possibly through lattice QCD as an intermediary. An alternative 
possibility will also be considered in this work. If the 'fundamental' theory at higher- 
energy is known and calculable, then at a momentum scale where the theories meet the 
parameters of the effective theory can be determined by matching onto the fundamental 
theory. This is sometimes done because certain calculations are more easily performed 
in the effective theory, either because of additional approximate symmetry in the low- 
energy limit or because a non-relativistic framework may be used. For instance this sort 
of matching calculation has been performed and applied with success in non-relativistic 
QED (NRQED) and non-relativistic QCD (NRQCD) 0. 

In this thesis we consider two applications of effective field theory to exploit the 
SU(3)l x SU(3)r chiral symmetry of QCD. In Chapter |2| we cover the theoretical frame- 
work upon which the effective field theories will be built. We discuss the symmetries 
of QCD which constrain the effective Lagrangian, the principles for constructing an effec- 
tive Lagrangian for the hadron degrees of freedom, and the power counting schemes that 
apply to the sectors of the theory with only light fields, one heavy field (static case), or 
more than one heavy field (non-relativistic case) . In Chapter [3] we present an effective field 



theory calculation of the self-energy shift of spin-| decuplet baryons in nuclear matter. 
With the exception of an expanded discussion of the A-baryon self-energy, the material 



in this chapter has already been published ||11|| . In Chapter ^ we present new material on 
an effective field theory for low-energy interactions of pions with kaons or an eta meson 
as an alternative to the standard SU(3) chiral perturbation theory. The objective is to 
achieve better convergence for very low energies by treating the kaon in a non-relativistic 
framework. To determine the parameters of the low-energy theory we perform a matching 
calculation from the heavy kaon/eta theory onto SU(3) chiral perturbation theory in the 
spirit of NRQED or NRQCD. Finally, we consider KK scattering in the heavy kaon 
formalism and discuss the utility of the matching calculation. 
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Chapter 2 

Theoretical Background 



In this chapter chiral perturbation theory (%PT) is introduced as the foundation for 
describing interactions involving the light pseudoscalar meson octet. As an effective field 
theory, xPT gets predictive power from the symmetries of the underlying theory (QCD) 
and a consistent scheme for counting powers of 'small' momenta. Section |2.1| reviews 
the relevant symmetries of the QCD Lagrangian. In sections |2.2| and |2.3| we outline the 



formulation of SU(3) xPT and its extension to include heavy fields. The momentum power 



counting for diagrams with one or two heavy particles is briefly discussed in section 2.4. 



2.1 Symmetries of QCD 

An appropriate starting point is the Lagrangian of QCD, 

£qcd = ~\G A ^G A ^ + q L ipq L + q R Hpq R - q R M q q L - q L M q q R + ■■■ (2.1) 

in which the ellipsis denotes gauge-fixing and ghost terms, renormalization counterterms, 
and the 8 term. As written £q CD is invariant under the Poincare group, SU(3)c gauge 
transformations, and charge conjugation (C). In addition the coefficient of the 9 term is 
known to be small so we neglect it throughout; in this approximation £q CD is also 
invariant under parity (P) and time reversal (T). 

In the limit that N of the quark masses vanish, the chiral limit, £q CD acquires addi- 
tional symmetries under independent U(iV) rotations of the left- and right-handed quark 
fields; 

1l q R -> Rq w (2.2) 



where L and R are unitary matrices restricted to acting in the massless sector. The 
U(l)y subgroup of this symmetry corresponds to the conservation of quark number in 
the massless flavors. A second subgroup, U(1)a, is broken in the quantum theory by 
the axial-vector anomaly. The symmetries of £q CD relevant for our purposes are the 
remaining SU(./V)l x SU(./V)r chiral symmetry and the discrete symmetries C, P, and T. 

The quark masses appearing in £q CD are non-zero and explicitly break the chiral 
symmetry. Even in the limit of massless quarks the vacuum expectation value of the 
quark bilinear 

(0|M ra |0) = -A5 afo (2.3) 

spontaneously breaks the chiral symmetry SU(A^)l x SU(./V)r down to the vector sub- 
group SU(-/V) V . The scale A x ~ 1 GeV, associated with spontaneous chiral symmetry 
breaking, determines the relative importance of the quark masses in breaking the chiral 
symmetry [pL3| . Quark masses m q which are much less than A x , specifically m u ^ < 10 MeV 



and 75 MeV < m s < 170 MeV |]14 |, can be treated as perturbations about the chiral limit 
m q — > by expanding in powers of m q /A x . For the rest of this chapter we identify the 
chiral symmetry group as G x = SU(3)l x SU(3)r and drop any explicit reference to the 
heavy quark flavors. 

Invariance under C, P, and G x is imposed on the Lagrangian of the effective field 
theory; invariance under T follows automatically from the CPT theorem. 



2.2 Chiral Perturbation Theory 

The spontaneous breaking of chiral symmetry to the SU(3)v subgroup implies the exis- 
tence of eight massless Goldstone scalar fields. However, because the chiral symmetry 
is explicitly broken by the quark masses, these fields acquire (small) finite masses and 
are commonly referred to as pseudo- Goldstone bosons. The pseudo-Goldstone bosons of 
spontaneously broken chiral symmetry are identified as the light pseudoscalar meson octet 
of pions (7r + , 7r~, and 7r°), kaons (K + , K~, K°, and K°), and the eta (rj). Chiral perturba- 
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tion theory is the effective field theory for describing the interactions of this meson octet 
at energies much lower than the chiral symmetry breaking scale, A x . For reviews of this 



subject see references [J5[ (Lf|. 

The pseudo-Goldstone bosons are represented by a 3 x 3 special unitary matrix, 
U = e**/ Fo , with 



( ^ 



7T 



7l n 



\ 



TV 

K~ 



-n"/V2 + r]/VQ K° 

-2 V /V6 J 



K° 



(2.4) 



and Fq ~ F n ~ 93 MeV. The pion decay constant F n is determined from 



{0\J»\*-{p))=iy/2F vJ t 1 , 



(2.5) 



where = wy^~f 5 d is an octet axial-vector current associated with chiral symmetry. Un- 



der the symmetries of section 2.1 the field U transforms as 



(2.6) 



We introduce the quark mass matrix A4 q through the field x — 2B Ai q , where the con- 
stant Bq is related to the vacuum expectation value in equation ( |2.3j ) and is approximately 
Bq ~ 2 A / Fq ~ 1300 MeV |16|] . For the purpose of constructing effective Lagrangians \ is 
assumed to transform in a way that preserves the symmetries of the QCD Lagrangian; 

X^RXL\ X ~^x\ X^X T - (2.7) 

In treating AA q we neglect the quark mass difference m u — and replace both m u and ma 
with the average rh = \{m u + m^). Because the corresponding SU(2) subgroup of SU(3)v 
is isospin symmetry, this approximation ignores isospin violation in the strong interaction. 
Calculations in x?T are organized as expansions in powers of m q and Q, the charac- 
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teristic momentum scale of the interaction. The terms in the Lagrangian of (£ xPT ) 
are grouped by the number of powers of m q and Q each contributes to diagrams when it 
is present. For a suitable choice of regulator and subtraction scheme, such as dimensional 
regularization with modified minimal subtraction (MS), each derivative contributes one 
power of Q. In the expansions of £ xPT (in powers of Q and m q respectively), the first 
terms are given by 

C xPT = *fTr[dUdW] + ^Tr[ X f/t + ^t] + . . . (2 . 8 ) 

The first term yields a canonically normalized kinetic term for the pseudo-Goldstone 
bosons <p a . The role of the arbitrary coefficient for the second term is played by the 
empirically-determined constant Bq. The second term gives the leading contribution to 
the pseudoscalar masses, 

ml = 2B m + ---, (2.9) 

m 2 K = B (rh + m s )-\ , 

m 2 = |5 (m + 2m s )H , 

represented collectively as m|. 

The result for the pseudoscalar masses indicates that the dual expansion in m q and Q 
is unnecessary. Since equation ( |2.9|) shows m\ ~ B Q m q ~ x an d we work in the relativistic 
regime where Q 2 ~ m%, one factor of m q (i.e., x) counts as a contribution of order Q 2 . 
Thus, the dual expansion is replaced with a scheme counting powers of Q only. The 
expansion of the Lagrangian is written £ xPT = C 2 + C i + ■ ■ ■ where C d contains all terms 
of order Q d and £ 2 is given by equation ( |2.8| ). The conventional parameterization of £ 4 
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1 


T r ( rn Will 3 


SO 111 Ce 


1 


0.4 ±0.3 


-K"e4, 7T7T — ► 7T7T 


2 


1.35 ±0.3 


-K"e4, 7T7T — ► 7T7T 


3 


-3.5 ± 1.1 


-K"e4, 7T7T — ► 7T7T 


4 


-0.3 ±0.5 


Zweig rule (large N c ) 


5 


1.4 ±0.5 


ratio Fjf : 


6 


-0.2 ±0.3 


Zweig rule (large N c ) 


7 


-0.4 ±0.2 


Gell-Mann-Okubo, with L r 5 , L r 8 


8 


0.9 ±0.3 


ratio (m s — rh) : (m^ — m u ), 









Table 2.1: Phenomeno logical values of the coefficients L%(/j), renormalized in MS 
at u, = m p , taken from J. Bijnens et al. |21] (also see reference |2(J for discussion). 



is 



£ 4 = L l Tr[d li Ud ft lf i ] + L 2 Tt[d fl Ud v U i ]Tr[d tt U& , U^] 

± L 3 Tr[d fl Ud"U^d u Ud u U^ + L A r R:[d ft Ud ft Ul]Tr\xU i + U X ] ] 
+ L 5 Tr[d^UdW( X Ui + C/X f )] + L 6 Tt[ x U^ + U X ] f 

+ L 7 Tt[ x U^ - urff + witixtfxrf + u x ^u x ^] 



(2.10) 



as first worked out by Gasser and Leutwyler |17|, [18|, [H], ^0[| . Table ^4] presents phe- 
nomeno logical values of the coefficients Lj renormalized in MS. In Appendix |A| we present 
Mathematica routines to symbolically expand £ x pt i n terms of the meson fields to the 
order necessary for this work. 

The power counting rules for Feynman diagrams, determined by Weinberg ||, establish 
the relative importance of diagrams in an arbitrary process. A diagram with Nl loops, 
N n meson propagators, and constructed from Nj vertices derived from C d will contribute 



at order Q where 
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D = 2 + 2N L + J2 N d( d ~ 2 ) (2- 11 ) 

d 

= AN L -2N w + Y,dN d . 

d 

The leading-order contribution is given by the sum of all tree diagrams constructed ex- 
clusively from operators in C 2 . At next to leading order one must include all one- loop 
diagrams constructed exclusively from operators in C 2 and all tree diagrams with one 
vertex from £ 4 and any number of vertices from £ 2 . Diagrams with more loops and more 
powers of Q 2 at the vertices must be included at higher order. 

On dimensional grounds, the powers of Q 2 which suppress higher-order contributions 
must be accompanied by a mass scale squared (A 2 ) in the denominator. For powers 
of Q 2 generated by a loop integral, the factor which appears is A 2 ~ (47r.Fo) 2 . When a 
factor of Q d arises from an operator in £ d , the compensating powers of A are implicit 
in the dimension of the coefficient of the operator. We can represent the typical size 
of coefficients appearing in C d as Cd = aA 4 ~ d , where a is a dimensionless constant of 
'natural' size (discussed below). Since each successive order in the expansion of £ xPT 
serves to approximate physics at short distance better than the preceding order, the 
scale relating coefficients in different orders of C xPT is characteristic of the short-distance 
physics approximated. Thus, A represents scales like the mass of the p meson (m p ) and 
the chiral symmetry breaking scale (A x ). We treat all of the scales which suppress powers 
of Q (e.g., 47rF , A x , m p , . . . ) equally in terms of the power counting and refer to the 
common scale as A = A x ~ AttFq. 

The convergence of calculations in depends on the assumption of naturalness, 
that the coefficients appearing in C xPT are not much larger than their natural size. An 
anomalously large coefficient, by a factor of order (A x /Q) 2 , would indicate that the cor- 
responding operator violates the power counting in equation ( |2.11| ) and the systematic 



expansion breaks down. A priori, there is no reason to expect such large coefficients 
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and the occurrence of such an anomaly would indicate either the omission of a relevant 
degree of freedom or fine-tuning in the underlying theory. We get an extremely simple 
estimate of the natural size of the dimensionless parameter a by recognizing C2 ~ Fq or 
a ~ (F /A x ) 2 ~ (47r)~ 2 . For instance, we expect the coefficients in £ 4 to be of about the 
natural scale, Lj ~ C4 = a. Manohar and Georgi |E| give a more rigorous derivation of 
the same estimate based on naive dimensional analysis. The phenomenological values 



presented in Table |2.1| are slightly smaller than this estimate. 



2.3 Coupling to Matter Fields 

The so-called matter fields are strongly interacting particles with masses which do not 
vanish in the chiral limit. Further, in the applications we consider, the masses of the 
matter fields are large enough that the fields can be treated in a non-relativistic formalism. 
The matter fields do not form representations of the full chiral symmetry group G x , 
but instead form irreducible representations under the approximate SU(3)v symmetry of 
the vacuum. In this section we describe a prescription for including the baryon octet 
fields (the N, S, and H isomultiplets and the A) and decuplet fields (the A, £*, and 
S* isomultiplets and the Q~) in the effective Lagrangian. For a review of this subject see 
reference [f[3[]. 

Under a G x transformation (L, R), fields in the fundamental representation of SU(3)v 
transform as ip — > Hip where H e SU(3) is a function of L, R, and the pseudo-Goldstone 
bosons U. The pure SU(3)v transformations correspond to taking H = L = R. The form 
of H is determined by specifying the transformation property of u = \/U = e % ®l 2Fo as 

u RuH ] = Hutf (2.12) 



which, when solved for H, gives 



h = y/mw lVw = VlwrY rVu. 



(2.13) 
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(The field u is an alternative representation for the pseudo-Goldstone bosons, correspond- 
ing to a different choice for the spontaneously broken generators of G x in the CCWZ pre- 
scription [p], [2^, 24fl .) Because the matrix if is a function of the pseudo-Goldstone boson 
field U(x), H also implicitly depends on the space-time coordinate x^. To compensate 
for the x^-dependence of H, we introduce a vector field 



(2.14) 



which transforms under G x as V — > HV^H^ — (d H)W and construct the covariant deriva- 
tive D ijj = d^ip + V J) such that DJ> -> H(D ip). 

Armed with the SU(3)v formalism for the fundamental representation, incorporating 
the spin-i octet and spin-| decuplet fields is relatively straight forward. The octet baryons 
are encoded as a 3 x 3 matrix of Dirac fields 



B 



f E%/2 + A/V6 £+ 

S- -T,°/y/2 + A/VE 



V 



WO 



V 



n 



\ 



-2A/V6 J 



(2.15) 



and transform in the adjoint representation of SU(3)v, i-e., 



B 



hbhK 



DB = dB + [V^B]. 



(2.16) 
(2.17) 



The spin-| decuplet baryons form a fully-symmetric rank-3 tensor under SU(3)v charac- 
terized by 



-\Q,bc &X j_j-LLLL J^ UfJ ^^ CU r J^ ]U 



T aa' TT bb' TTCc'ma'b'd 



(DT) 



abc 



dm&bc i t r&a> rn& be . T rbb m&b c . T t-cc m® 
T +V T +V T +V T 
ii fi fi fi 



(2.18) 
(2.19) 



The components of T are Rarita-Schwinger fields subject to the auxiliary constraint 



7^T = and are identified as follows: 
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T m = A++, T 113 = Z*+/V3, T 133 = E*°/V3, (2.20) 

T 112 = A+/V3, T 123 = E*°/v / 6, T 233 = E*~/y/3, 

T 122 = A°/V3, T 223 = E*"/V3, 

T 222 = A-, T 333 = Q-. 

The minimal-coupling Lagrangian for the -B and T fields is 

£ = B(iip-m B )B-T fi ( y iip-m T )T^ (2.21) 

with implied pair-wise summation of the dangling chiral indices. 

The effective Lagrangian for matter fields will also explicitly include U and \ which 
transform under G x as O — > ROL^ . For convenience, we define fields O = u'Ou* such 
that O — > HOW under G x transformations and coupling B and T to the new fields O 
becomes transparent. The result for U is trivial, U = u'Uw = 1; for d U we get the more 
interesting result 

\ = WW* 1 = T^u (2.22) 
= '3 u - ud v)) , 

where the factor of | is included so is hermitian. To replace \ (and x^) we choose the 
(anti-)hermitian combinations 

X± = u X u =t ux'u. (2.23) 

As an example the order-Q 2 Lagrangian of x?T, equation ( |2.8D , can be rewritten in terms 
of and X ± as C 2 = F 2 TrL4^] + fF 2 Tr[ X+ ]. 

We briefly mention two relations which help to identify a minimal set of terms which 
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are needed at higher order, see Fearing and Scherer 
The first is the chain rule in the form 



231 



for a more complete account .Q 



■ft 



^ = 50-^,0} (2.24) 



which allows the replacement of multiple derivatives of U with factors of and co- 
variant derivatives of A , e.g., i$d d u Uv) = —2iD A v — 2{A^A u }. Two consequences of 
equation ( |2.24|) are 1) D A v is symmetric in (ji, v) and 2) despite the dependence of %± 



on u(x), covariant derivatives of x± are unnecessary since D x± — i{A , Xt\- The secon d 
simplifying relation gives the field strength associated with in terms of the field A^ 

G (v) = d V -d V + [V ,V] = [A ,A ]. (2.25) 

Because antisymmetric covariant derivatives result in factors of the field strength, like 
{D^D V — D y D^)ip = G^Jij) in the fundamental representation, we can treat all multiple 
covariant derivatives of any field as implicitly symmetric without loss of generality. To 
summarize, the G x -invariant effective Lagrangian for the baryon fields B and T is written 
in terms of the building blocks B, T, A , x±i an d fully-symmetric covariant derivatives 
of them. The field strength G^ V J and covariant derivatives of x± may be omitted in favor 
of terms involving more factors of A . 



2.4 Non-Relativistic Power Counting 



The appearance of baryon masses in equation ( p. 21 ) wrecks the power counting of xPT 



Derivatives of heavy fields would contribute factors of 'hard' momenta, where p > ms, 
and loop diagrams could result in explicit factors of mg or in the numerator. These 
diagrams are not suppressed relative to 'lower order' diagrams since mg ~ itlt ~ A x . 

We use a non-relativistic approach adapted from the heavy baryon chiral perturbation 

^^Whcn using the external field method to impliment local chiral symmetry, these relations are modified 
by additional terms involving external (i.e., non- propagating) gauge fields. 
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theory (HB%PT) developed by Jenkins and Manohar |26], |27|, [29| to address this 
problem. We start with the velocity-dependent baryon fields of HB^PT then specialize 
to working in the rest frame of the baryons. 

The momentum of a heavy baryon field is decomposed as = msv^ + k^ in HB;\BT, 
where the residual momentum k^ is assumed small and reflects how much the baryon is 
off mass-shell. The velocity-dependent baryon fields are defined using v M as 

B v (x) = ^±t e imBV - x B(x), (2.26) 
T v {x) = ^±t e imBV - x T{x), 

where the octet mass is used in the definition of the decuplet T v to avoid ^-dependent 
phases in the Lagrangian coupling B v to T v . The factor of |(1 + f) projects out the particle 
components of the Dirac spinors. The anti-particle components are implicitly integrated 
out of the theory and the effects of virtual baryon loops are absorbed into terms of the 
effective Lagrangian suppressed by powers of l/m#. 

This representation of the heavy baryon fields permits a sensible power counting 
scheme. Derivatives of the velocity-dependent fields give factors of the small residual 
momenta k^ in place of the hard momenta p M . The minimal-coupling Lagrangian corre- 



sponding to equation ( |2.21|) becomes 



C v = B v (iv ■ D)B V - T v ^(iv ■ D)T/ + AmT^T/ + • • • , (2.27) 

where Am = my — ttib — 300 MeV is considered of order Q and the ellipsis denotes 
higher-order terms, such as —B v D 2 B v /2rriB, induced by the integration over anti-particle 
degrees of freedom. By removing mg and tut from the baryon propagators, loop diagrams 
will not introduce positive powers of the masses in diagrams (except in a case considered 
toward the end of the section). 

The Dirac structure of the fields can be eliminated in favor of two-component spinors 
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because we have integrated out the anti-particle components. We choose to explicitly work 
in the frame v = (1,0,0,0) and drop the subscript v on heavy fields. This choice of frame 
simplifies the spinor-related notation; for instance, the auxiliary condition on Rarita- 
Schwinger fields, 7 = 0, implies constraints which reduce to T° = and a • T — 0. In 
this frame the resulting non-relativistic framework is equivalent to a Lagrangian formu- 
lation of the time-ordered approach discussed by Weinberg [HJ . 

The power counting for diagrams with a single heavy field generalizes equation ( |2.11| ). 
Since uib ~ rnr ~ A x we do not need to distinguish between corrections suppressed 
by Q/rriB and Q/A x . As before, each meson propagator counts as Q~ 2 and each loop 
integration gives a factor of Q A . Heavy field propagators each contribute Q~ x as seen 
from equation (|2.27|) . We let Ni represent the number of vertices in a diagram which con- 



tribute di factors of Q and contain rij heavy fields. A diagram in the single heavy-particle 
sector with N L loops, N n meson propagators, and Nj baryon propagators contributes at 
order Q D where Q 

D = 4N L -2N W - N T + J2 d i N i ( 2 - 28 ) 

i 

= l + 2iV L + ^iV i ( C / l + in i -2). 

i 

The primary difference from equation ( |2.11| ) is that each order is suppressed by Q/A x 
rather than (Q/A x ) 2 relative to the preceeding order. 

For diagrams with two or more heavy particles, the power counting is complicated 
by infrared divergences in some loop integrals ||31|| . Consider a heavy-particle bubble 
diagram; for octet fields the loop integral takes the form 

Y~Y — f-^- i % — = — ! dd ' lq (2 29) 

J(^) d E + q +ie E-q + ie 2E J (2n)^ 7 K ' } 

where the incoming particles have energies E and momenta ±p. The divergence arises 
as E — > because the q contour is pinched between the poles of the static propagators. By 
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resumming the kinetic energy operator, e.g., B^V 2 B/2mB, into the heavy field propagator, 
the infrared divergences are removed. In dimensional regularization the result for the 
modified loop integral becomes 



d d q i i 

{2ii) d _|_ g Q _ \q\ 2 /2rriB + it E — q — \ — q\ 2 /2rriB + ie 

d d ~ 1 q i 
(2vr) d - 1 2E - q 2 /m B + 2ie 

rf-3 

irriB f—2m B E — ie'\ 2 1 _, /3 _ d s mB\/^ m BE 



(2.30) 



47T V 47T y/ 2 ' 47T 

which is well-defined in the infrared. Setting E = p 2 /2rriB + ■ ■ ■, the leading behavior of 
the loop integral is ~msQ, not ~Q 2 as expected from equation Q2.28|) . 



By including the kinetic energy in the heavy field propagator we are treating the 
operator V 2 /2m# on an equal footing with d Q . For consistency, powers of the integra- 
tion variable q Q should be counted as q ~ Q 2 /ms in loops leading to nearly-infrared- 
divergent behavior. Heavy field propagators count as tub/Q 2 , loop integrals now give 
dq Q d 3 q ~ Q 5 /m B , and meson propagators contribute 1/Q 2 as before. Thus the loop in- 
tegral above counts as ~(Q 5 /^s)(^s/Q 2 ) 2 = ttlbQ which reconciles the power counting 
scheme with the result of equation (|2.30|) . A general Feynman diagram with heavy par- 
ticles in the initial and final states may contribute at order Q D where 

D = 5N L - 2N W -2N! + Y1 d * N * ( 2 - 31 ) 



with the same notation as used in equation ( 2.28| ). If the loop integrals in a diagram are not 



infrared divergent, as with crossed pion exchange, the actual contribution will be of higher 
order consistent with equation Q2.28p . The systematics of the power counting for diagrams 
with two or more heavy particles is developed in references 0, [31], [32], |34], [35], [36]. 



The presence of a bound state or resonance near threshold can also complicate the power 



counting |}7, q^, |39|. This complication arises in nucleon-nucleon scattering and the 
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power counting was studied in that context by several authors, in particular Kaplan, 
Savage, and Wise 0, |41] . In this work we do not consider systems requiring this special 
treatment and refer the reader to the literature for a discussion of the relevant power 



counting schemes H2L El IP 
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Chapter 3 

Decuplet Self-Energy in Nuclear Matter 

Strong interaction effects shift the self-energy of hadrons in nuclear matter from the free- 
space values. This effect for the spin-| octet baryons was studied in an effective field 
theory framework by Savage and Wise [45]. Using the formalism reviewed in Chapter ^ 
we calculate the self-energy shifts of the spin-| decuplet baryons in nuclear matter at 
leading order. The self-energy shifts of the decuplet baryons, particularly of the A iso- 
multiplet, are relevant in studies of meson- nucleus scattering |46| , |47fl and of stellar and 
neutron star matter |48| , [4~9| . For the A isomultiplet the self-energy shifts have also been 
examined in various phenomenological models |5D|, ^TJ and in QCD sum rules [|52|| . 
The work described here differs from the earlier approaches by extending the calculation 
to multiplets of chiral SU(3)v and by including contact diagrams necessary for a consis- 
tent and systematic momentum expansion. This chapter makes a minor correction to a 



prior publication JTTJ] and includes a more detailed discussion of the self-energy of the 
A isomultiplet in nuclear matter near saturation density. 

The first section discusses how effects of nuclear matter are described in the effective 



field theory. In section |3.2| we determine what Feynman diagrams contribute to the self- 



energy shifts at leading order and construct the relevant effective Lagrangian. Section 3.3 



presents the main results of the calculation. Finally, a discussion and interpretation of 



the results is contained in section 3.4 



3.1 Effects of Nuclear Matter 



In nuclear matter the propagation of decuplet baryons is effected by interactions with 
the background medium. In particular the self-energy S, i.e., the location of the k pole 
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in the decuplet two-point function, is shifted relative to free space. As an additional 
consideration, the background medium breaks Lorentz boost invariance by specifying a 
unique frame, the zero-momentum frame of the nuclear matter. The self-energy in free 
space of a decuplet baryon with momentum k — \k\ is 

£ vac (> 2 ) = \J m T + k 2 - im T T vac - m B (3.1) 

= ( Am + ^ + ''H r ™( 1 -S + ''') +0(r " J ' 

where T vac is the free-space decay rate and Am = m T — m B . Throughout this chapter, 
the octet mass mg has been implicitly subtracted whenever we refer to the decuplet 
self-energy, as discussed in section ^j4]. Equation (|3~T|) contains only two decuplet param- 
eters; the full fc 2 -dependence of the self-energy is determined by Lorentz invariance. The 
corresponding expression for a decuplet baryon in nuclear matter can be written 

£ nm (k 2 ) = (Am* + ak 2 + ■■■)- §r*(l + (3k 2 + ■ ■ ■) (3.2) 

in which the /c 2 -dependence (e.g., a, (3, . . . ) cannot be determined from symmetry argu- 
ments alone. 

Although invariance under Lorentz boosts is lost, the remaining rotational symmetry 
constrains the spin-dependence of the self-energy £ nm . Because the nuclear medium is 
rotationally invariant, the only preferred spacial directions are along the decuplet baryon 
three-momentum k and spin S. As a function of k and S, the self-energy depends on 
only the combinations k 2 and k ■ S. (S 2 = ^ is trivial.) Further, parity invariance of 
the strong interaction dictates that the self-energy depends on even powers of k, which 
means replacing k ■ S with (k ■ S ) 2 = k 2 h 2 in terms of the baryon helicity h. Consequently, 
the self-energy in nuclear matter is diagonal in the baryon helicity states and takes the 
values ^ 2) (fc 2 ) for h = ±| and ^ 21 (/; 2 ) for h = ±|. In the limit of vanishing momentum, 
the rotational symmetry is elevated to full SU(2) invariance, and the self-energy must 
be independent of the decuplet spin projection along any direction. In terms of the 
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parameters in equation ( |3.2| ), Am* and T* are helicity- independent, while the coefficients 
of powers of k 2 (e.g., a, /3, . . . ) depend on \h\. 

The lowest order in a density expansion for nuclear matter is a Fermi gas of non- 
interacting protons and neutrons with Fermi momenta p ( p and p { p respectively. In this 
framework, the characteristic momenta relevant in the chiral derivative expansion are p£ , 
Pp\ and k the decuplet baryon momentum. Since the density of a degenerate Fermi gas 
is given by dp = p F /3^ 2 , the density expansion for nuclear matter is consistent with the 
chiral derivative expansion. The static nucleon propagator in nuclear matter with Fermi 



momentum pp is [45] 



Snm{q ,q) = + (3.3) 

at lowest order in the nuclear density. The modified nucleon propagator reflects the 
presence of the background medium through two effects, nucleon states inaccessible due 
to Pauli-blocking and nucleon- hole intermediate states allowed for \q\ < pp. 

As the location of the k pole in the exact decuplet two-point function, the self-energy 
in nuclear matter is given by the solution of 

S nm (k 2 ) - Am - E nm (£ nm (fc 2 ), k ) = 0, (3.4) 

where S nm is the proper self-energy for nuclear matter, i.e., —iT> nm is the sum of connected 
one-particle-irreducible diagrams in the two-point function. What we calculate is the 
self-energy shift, 58 = £ nm — £ vac , obtained from equation ( |3.4| ) by expanding the proper 
self-energy for nuclear matter about the free-space pole £ vac , 



5£(k 2 ) = E nm (£ vac (/ C 2 ),^)-E va ,(£ vac (/ C 2 )^) (3.5) 
+ {^Jk 2 )-S vac (k 2 )}^UK,k] 



+ ■ ■ 



and in turn expanding £ yac (k 2 ) in powers of Q as shown in equation (^.1[). Note that 
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(a) 




Figure 3.1: Feynman diagrams for the self-energy shifts at leading order of decuplet 
baryons in nuclear matter, (a) meson-nucleon loop diagrams and (b) contact diagrams. 
Double lines represent decuplet baryons, single lines represent nucleons, and dashed lines 
represent pseudo-Goldstone bosons. 

only Feynman diagrams with an internal nucleon propagator contribute to the difference 
on the first line. Because the chiral expansion of the proper self-energy S nm starts at 
order Q 2 , successive terms in equation ( |3.5| ) only contribute at higher order. (A very 
important exception to the last point is discussed in detail in subsection p.4.2 .) The real 



part of the self-energy shift 5E (h) (k 2 ) modifies the decuplet baryon energy-momentum 
dispersion relation from the free-space form. The change in the decuplet decay rate is 
given by 5T (h) (k 2 ) = — 2 lm[5£ (h) (k 2 )} at leading order. 



3.2 The Effective Lagrangian 

At leading order in the chiral expansion the self-energy shifts 58 coincide with the 
difference £ nm (£ vac , k ) — £ vac (£ vac , k ). For diagrams contributing to the proper self- 
energies E nm and E vac , the power counting is given by equation ( j2.28|) for the single 



heavy-particle sector. Because a nucleon propagator is required for a non-zero differ- 
ence, the leading contribution arises from one-loop diagrams in which the vertices satisfy 
d + |n — 2 = and contributes at order Q 3 . Topologically, the one- loop diagrams are con- 
structed from either two three-leg vertices or a single four-leg vertex; the two possibilities 
with a nucleon propagator are shown in Figure 

The vertices of meson-nucleon loop diagrams are derived from operators coupling the 
fields TB^A^ in the order-Q Lagrangian. The most general Lagrangian of that form, 
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invariant under G x and parity, is p7| , ^ 



£ a = - C e abc (A" 4 ■ [B^ be f cde ] + A da ■ [f t cde 5 efe ]) , (3.6) 
where square brackets denote summation on spinor indices. The value of the coefficient, 



\C\ ~ 1.53, is empirically determined from T — > Bn decays The meson-nucleon loop 
diagrams contribute to the self-energy shifts of the A and S* isomultiplets only; the 
Lagrangian £ a does not couple 5* or Q~ baryons to a nucleon and single pseudo-Goldstone 
boson, a consequence of strangeness conservation. 

For contact diagrams the vertex contributes d = powers of Q, so the relevant La- 
grangian contains only simple products of field operators TB{TB)^ . To construct the 
four-baryon operators we start with the spin and chiral structures of the product TB. 
The operator product decomposes under rotational SU(2) as | © | = 2 © 1 and under 
chiral SU(3) V as 10 © 8 = 35 © 27 © 10 © 8. By coupling products TB and (T5) t to 
form chiral and rotational singlets, we find eight linearly-independent four-baryon oper- 
ators which contribute to the self-energy shifts at leading order. We choose to write the 
effective Lagrangian which contains these operators asQ 



£ dl \rp]j abc rpjabc-i ^-\ed gde-i <^5 ^rp^jabc^k rpj abc-i ^g^ed^k gde-i (3 7) 



^2 ^r-p^j abc rpj abd^j^^ed j^ce-^ ^6 ^p,1[j abc^k rpj abd^p^jed^k j-^ce^ 



2F 2L JL J 2F 2 

d^ ^rp^j abc rpj abd^j^^ce p^ed-^ dj \^p]j abc ^k rpj abd-^ ^gj-ce^.fe^e<ij 
d^ rrriij abc rrij ade-i r j-.fbd j~.ce-! dg rmfjabc k mj ade-i r r^fbd k t-cc-i 

where {j,k) are vector indices, (a-e) are chiral indices, and square brackets indicate sums 
over spinor indices. Factors of F ~ 2 are included in C h to make the coefficients di dimen- 
sionless. 



1 In our prior publication we used a different convention for the pion decay constant, specifically 
/ = V2Fo ~ f 32 MeV. 
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The values of the eight coefficients dj have not yet been experimentally determined. 
However, because the baryon helicity is conserved in the self-energy diagrams, the terms 
in C h of the form [T^cr k T^][B^a k B] do not generate self-energy shifts and the results 
are independent of d^-dg. Values for the remaining four coefficients d\-d± are important 
for quantitative predictions. Unfortunately, knowledge of the values will likely have to 
wait until low-energy octet-decuplet scattering data become available. In section |3]4] we 
discuss constraints on the coefficients di from the further assumption of SU(6) spin-flavor 
symmetry. 

Two characteristics of the contact diagrams allow some predictions which are inde- 
pendent of the coefficients di. In a contact diagram there is no 'intermediate state' which 
corresponds to an allowed decay of the decuplet baryon so the di do not appear in the 
imaginary part of the self-energy shift. Also the contact diagrams are independent of the 
momentum on the external line so the c/j-dependence is restricted to the fc 2 -independent 
parts of the self-energy shifts. In terms of the parameters in equation Q3.2|) , at leading 
order only the helicity- independent quantity Am* depends on the coefficients d{. Quan- 
tities independent of the coefficients di, for which we present quantitative results, are the 
helicity-splitting of the self-energy shifts A^E = 5E (1/2) — 5E (3/2) and the decuplet decay 
rates in nuclear matter = r var + 5T (h) . 



3.3 Self-Energy Shift Results 

From the effective Lagrangians, equations ( |3.6| , |3T7| ), we calculate the self-energy shifts SE (h 
and ST (h) of the spin-| decuplet baryons in nuclear matter to leading order (~Q 3 ) in the 
chiral momentum expansion. For convenience we introduce the two 'threshold' mass scales 

H = JirriA -m N ) 2 -ml ~ 255 MeV, (3.8) 
\x = \ m 2 K — (m s * — m N ) 2 ~ 215 MeV, 
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and present the results in terms of the functions 



gi(k,p,m) 
Gk(k,p,m) 

2 

g3(k,p,m) 
Gz(k,p,m) 



—L^( m 2-j? + k *) f (m 2_ p 2 + A;2) 2 + 4m 2 k 2} 

lorrr/r I ) 

l — {3(m 2 -p 2 + k 2 ) 2 - 2k 2 (p 2 + 3k 2 -9m 2 )} 
2 -(p- kf 



(3.9) 



12m 2 k 2 



m 



H gi(k,p, m) In 

p 2 

-1 



{p + kf 



m 2 -p 2 + k 2 ) {{m 2 -p 2 + k 2 ) 2 - 12m 2 A; 2 | , (3.10) 



16m 3 A; 3 

I {S(m 2 -p 2 + k 2 ) 2 - 2k 2 (p 2 + 3k 2 + 15m 2 )} 

m 2 — (p — k) 2 



m 



H gz(k,p, m) In 

p 2 



(p + ky 



which are constructed such that Qh(k,p,m) — > as k — > + . 

For the negatively-charged member of each isomultiplet we find the real parts of the 
self-energy shifts are 



5E™ 



{(3d 1 )pf+ (3d, + 3d 2 )pP 3 } 



(3.11) 



+ 



48^r^ 3ln 



(fi-p { ; y ) 2 -k 2 

(n + pf) 2 -k 2 



+ ^f+Afi 2 p^ + ^ 2 ppg h (k,pp,A, 



5E™_ = ^^{(3d 1 + d 3 )p ( f+(3d 1 + 2d 2 + d 3 + d 4 )p ( f} (3.12) 



+ 



o 

C 2 



tf- p f 2 + tf 



— 77 < 2a arccos / 



5E, 



Q- 



{(3d, + 2d 3 )pf° + (3d, + d 2 + 2d 3 + d,)p { f] , 
1§ l^(3d 1 + 3<fe) {^ p)3 +^ )3 }. 



(3.13) 
(3.14) 
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The helicity splitting of the A" and £*" self-energies are graphed as functions of p ( p and k 
in Figure jO) . 

At leading order the imaginary part of the self-energies, or resonance widths, are 
shifted in nuclear matter for the A isomultiplet only. We find the A - resonance widths 
are 



r nm,A- = {\{l+9H{kJ;\^))Q{k+pf-^)Q{^-{pf-kf) (3.15) 

+ e( / *-j# ) -*) + (A;-^ ) - / ,)}r vac>A _ 

where T vacA _ = // 3 C 2 /127tFq is the A" resonance width in free space at leading order. 
Figure |3.3| presents the resonance widths as functions of pp and k. 



Results for the other baryons in the decuplet can be determined from the expressions 
for the negatively-charged members. For the member of each isomultiplet with the most 
positive charge (A ++ , E* + , and the self-energy shifts are obtained by exchanging pp 
and pp in equations ( |3.11| - |3~T3| , |3.15| ). The self-energy shifts of the remaining decuplet 



baryons are given by the following relations: 

5£ A+ = \{28£ A++ + 8£ A -), (3.16) 
5£ A o = l(6£ A++ + 25£ A -), (3.17) 
5£x*o = l(6£ E , + + 5£^-) . (3.18) 

Clearly, in neutron-proton symmetric nuclear matter the baryons within isomultiplets 
remain degenerate up to the splitting between helicity states. 



3.4 Discussion and Conclusion 
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Figure 3.2: Leading-order helicity splitting of the decuplet baryon self-energy in nuclear 
matter A (h) E = 5E {1/2) — 5E is/2) as a function of neutron Fermi momentum and baryon 
momentum k, (a) for the A" baryon and (b) for the £*~ baryon. 




Figure 3.3: Leading-order A' resonance widths in nuclear matter as functions of neu- 
tron Fermi momentum pp } and baryon momentum k, (a) for helicities ±| and (b) for 
helicities ±|. 
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(b) 



Figure 3.4: (a) Infrared divergent baryon-meson box diagram and (b) overlapping meson- 
loop diagram obtained from the box diagram by contracting nucleon lines 2 and 3. The 
two- loop diagram in (b) avoids the infrared divergence of (a). 

3.4.1 Corrections at Higher-Order 

Contrary to the discussion in our prior work, the leading corrections to our results are 
not from infrared divergent two-loop diagrams. We previously thought diagrams with 
overlapping meson loops would be afflicted with the infrared divergence that appears in 



baryon-meson box diagrams, see Figure |3.4| . That the two-loop diagrams avoid the di 



vergence is most easily seen in terms of time-ordered perturbation theory as discussed by 
Weinberg [^, |3~1J| , where field propagators and integrals over four- momenta are replaced 
by intermediate-state energies in the denominator and integrals over three-momenta. 
In the time-ordered approach, infrared divergences arise from intermediate states with 
'small' energies of order Q 2 /rriB', however, in loop diagrams contributing to the proper 
self-energy H(k Q ,k) all intermediate states contain mesons with energies on the char- 
acteristic scale Q. Consequently the two- loop diagram in Figure |3]4](b) is governed by 
equation ( |2.28| ) and is suppressed relative to the leading-order result by (Q/A x ) 2 and not 
by m B Q/A 2 x ~ Q/A x . 

The leading-order results (~Q 3 ) presented in section |3^ are determined from the 
difference E nm (Am, k ) — £ vac (Am, k ). A priori, we would expect corrections at order Q 4 
from three sources: 1) the chiral expansion of the proper self-energies S nm and S vac 
based on the power counting of equation (|2.28| ), 2) the expansion of S nm about £ vac in 
equation (|3.5D, and 3) the Q-expansion of £ vac as £ vac = Am + k 2 /2niT + 0(Q 3 ). In the 
following two paragraphs we argue that the corrections from the first two expansions do 
not occur until order Q 5 . The only corrections of order Q 4 arise from the third expansion 
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and are included by making the replacements 



/Am\ k 2 

"^" + (— )^r T (319) 



[l > [L 



Am\ k 



,2 



jj J 2rriT 



in the leading-order results for the A and E* isomultiplets. The results for the S* isomul- 
tiplet and the Q~ baryon do not receive corrections at order Q 4 . 

The chiral expansion of the proper self-energies S nm and S vac would contribute at 
order Q through one-loop diagrams obtained from Figure ^Tl] by replacing a single vertex 
with one which satisfies d + ~n = 3. The only relevant operators, invariant under rota- 
tions, P, and G x , are constructed from products of the form (id )TB*A for the meson- 
nucleon loop diagrams or (id Q )TB(TBy for the contact dia grams. Partial integration and 
baryon field redefinitions can be used together to eliminate both classes of operators in 
favor of terms higher-order in the expansion of the effective Lagrangian. (Such field re- 
definitions are reviewed in Chapter [|.) Consequently, there are no corrections at order Q 4 
from the chiral expansion. Order-Q 5 corrections arise from one-loop diagrams with a 
single vertex which satisfies d + ~n — 4 and two- loop diagrams, such as in Figure ^4|(b). 

The primary corrections from the expansion of E nm about £ vac are suppressed by 
powers of Q contributed by 9S nm /9fc in equation (^.5[). Because the Q-expansion of S nm 
begins at order Q 2 , e.g., through an insertion of the kinetic energy operator, we would 
expect the corrections to be suppressed by Q, or to contribute at order Q 4 . However, all 
/c Q -dependence of S nm can be moved into loop diagrams which are order Q 3 or higher by 
suitable redefinitions of the decuplet baryon field. Thus, <9£ nm /<9fc ~ Q 2 and the primary 
correction contributed by the expansion of S nm about £ vac is order Q 5 . 

Throughout this calculation we have made the assumption that the power counting in 
equation (|2.28|) applies. However, large scattering lengths for nucleon-nucleon or nucleon- 



delta interactions will require a revised power counting schemer] as seen in references [HO 



2 We thank M. Wise for bringing this point to our attention. 



29 



4l|| . Using the revised power counting for the problem considered here is presently not 
tractable. How these effects may be taken into account remains to be seen; the work 
presented here is the best that can be done at the present time and uses the power counting 
of equation ( |2.28| ). We hope that the discussion of some technical details, particularly 
those presented in the next subsection, will be useful in further studies of the decuplet 
self-energy. 

3.4.2 A Self-Energy Near Saturation Density 

One important feature of the self-energy shifts of the A isomultiplet, not apparent in 
Figure pT2| (a) , is a logarithmic divergence when k = at a 'threshold' Fermi momentum 
(coincident with the discontinuity of the widths). The divergence arises due to degeneracy 
of the A baryon with a pion and a nucleon on the surface of the Fermi sea; for static 



baryons p* F = fi, which is just the solution of wa = + -J m\ + pp . To the extent 
that the divergence complicates the analysis of the self-energy shifts, it is unfortunate 
that p* F ~ 255 MeV is so close to the Fermi momentum associated with nuclear matter 



at saturation density, p^ at ^ ~ 262 MeV [j53j. In contrast, the E* isomultiplet self-energy 



shifts do not have a similar divergence because the X* — > NK decay is kinematically 
forbidden for any p { p and Pp) . 



In the A self-energy shift, equation (|3.11|), the divergence manifests itself as a can- 



cellation among quantities of order Q in the argument of the logarithm, i.e., when 

Ue' 1 -m\- pp) =fi -p F = 0. (3.20) 

V / £=Am 

The divergence is resolved by keeping higher-order terms in the argument of the logarithm, 
specifically the imaginary part in the expansion of S about Am. When the imaginary 
part is kept, the logarithmic divergence is replaced by a logarithmic enhancement of 
the self-energy shift and the result is order Q z \\iQ 2 . For Fermi momenta sufficiently 
below p F , the cancellation is not significant and we expect equation ( |3.11|) to provide 
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a good description of the A" self-energy shift. The effect of the cancellation becomes 
important when \i — p F ~ Q 3 , or approximately for p F > (0.7)p* F ~ 180 MeV. 

The logarithmically divergent term can be traced to the proper self-energy in nuclear 
matter E nm (£, 0). At the threshold Fermi momentum p* F , S nm is singular when £ = Am 
and the expansion in equation (|3.5| ) fails. To calculate the self-energy £ nm in the vicinity 
of p* F , we evaluate S nm at £ nm and expand £ vac (£ vac , 0) about £ nm (S vac is smooth in this 
region) , 



d — ^vac + ^nm(^nm5 0) ^vac(£ nm ,0) (3-21) 



(^nm ^vac) ni- ^vac(^O'O) 



'0 



+ 



k —f 

"'O nm 



As in equation ( |3.5| ), all the terms of order Q 3 appear on the first line and the second line 
may be neglected at leading order. 

In equation ( [3.2 1| ) we pay a price to determine the self-energy £ nm near the threshold 
Fermi momentum; £ nm is the solution of a transcendental equation and must be found 
numerically. Figure |3.5| presents a numerical calculation of the self-energy shift SE A _ com- 
pared with the divergent behavior of equation ( [3.11| ). Both curves are plotted as functions 



of p F with k = and we assume the coefficients of the contact terms are d\ = di = 0. Be- 
cause the contributions of the contact diagrams are effectively omitted, the results plotted 
in the figure are not model-independent in the sense of effective field theory. Until the 
values of the coefficients di and di are determined, there is little motivation for further 
numerical study of equation ( |3.21| ). 



Having resolved the logarithmic divergence in the real part of the A self-energy shift, 
next we consider the behavior of the resonance width r nm A _ near p* F . The basic features 



of the plots in Figure |3]3| are easily understood in terms of the nucleon momentum in 
A — > Ntc decay. When the momentum of the A~ baryon is k, then up to recoil correc- 
tions the nucleon and pion momenta in the final state are respectively k + ft and —jl, 
where equals the mass parameter fi. Naively, when \k \ + [A < p F all the nucleon final 
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Figure 3.5: Self-energy of the A" baryon at rest in nuclear matter, omitting contact 
diagrams (d\ — di — 0). The solid curve is the numerical solution of equation (|3.21|) ; the 
dashed curve is a plot of equation fl3.11Q . 



states are occupied by nucleons in the Fermi sea; the A" baryon becomes stable against 
A — > Nir decay and the width in nuclear matter vanishes. When pp < \\k\ — /i| all the 
nucleon states are available for decay and the width in nuclear matter is the same as the 
free-space width r vac . For values of pp between those two cases, the Fermi sea partially 
obscures the shell of nucleon momenta and < r„ m A _ < r,„„. When k — the width 

nm,/A vac 

transitions abruptly from r vac (for pp < p) to zero (for pp > p) creating the discontinu- 
ities in Figure j3~3] . 

The pitfall in the simple analysis of the previous paragraph is determining the momen- 
tum transfered to the intermediate-state nucleon p in terms of the free-space energy of 



the A baryon. More rigorous consideration of equation (3^) shows that the A resonance 
width is non-zero for any p F , although it may decrease sharply near p* F . Setting k = 0, 
£„ m is the solution of 



Am 



£ S nm (£) 
2Atx'< 



'iim v"^ J 

l-£ 2 )j \ f-e\ 



2 arctan 



[ Pf 



(3.22) 



+ f(S,p F ), 



which is analytically continued to £ > m n by assigning m^. an infinitesimal imaginary part, 



mi 



— ie. The last term, f(S,pp), is an analytic function of S below the threshold 
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for Nirn production from a virtual A baryon and does not influence our reasoning. The 
only real- valued solutions of equation (|3.22|) satisfy either £ < m n or p F > \J £ 2 — ra\ ; 
however, as pf is decreased with Am fixed, none of these solutions flow continuously to 
the physical, complex-valued solution for small pp. Starting with £ vac at pp = 0, as pf 
is increased the physical solution bypasses the region of real-valued solutions by passing 
onto a different sheet of the Riemann surface. 

The main result of this subsection is the logarithmic divergence of 6E A _ and the 

° ° nm,A 

discontinuity of r nm A _ , seen in section |D| can be removed by changing how the evalu- 
ation of S nm is handled. Quantitative predictions for £ nm near p* F , by numeric solution 
of equation ( |3.21| ) or ( |3.22| ), require values for the coefficients d\ and c?2- Independent 



of d\ and g^, however, we conclude that the resonance width of the A isomultiplet is 
non-zero for any pp] the decay channel A — > Ntt is never fully blocked by the presence of 
the nuclear medium. 

3.4.3 Conclusion 

To complete the description of the self-energy shifts of the spin-| decuplet baryons in 
nuclear matter, the coefficients dj of the contact terms must still be determined. In the 
absence of low-energy octet-decuplet scattering data, one way to estimate the values of 
the coefficients may be through appealing to a larger symmetry group, in this case the 
approximate SU(6) spin-flavor symmetry. Because the spin-| octet and spin-| decu- 
plet baryons form a single 56-dimensional representation of spin-flavor SU(6), the SU(6)- 
invariant Lagrangian determines the octet-decuplet coefficients di-d 8 , the Savage- Wise 
coefficients Ci-c 6 for octet-octet interactions |45[| , and 14 coefficients for decuplet-decuplet 



interactions. The SU(6)-invariant Lagrangian contains only two four-baryon contact terms 
of dimension six, with coefficients a and b as defined by Kaplan and Savage |54] . In terms 



33 

of a and b, the octet-decuplet coefficients in equation ( |3.7| ) are 



d 1 = 


2a + §6, 


d 5 = 




d 2 = 


9"' 


d 6 = 




d 3 = 




dj = 




C?4 = 




d 8 = 





(3.23) 



Unfortunately, neither a nor b is reliably known; therefore, for the phenomenologically- 
interesting A isomultiplet, we have exchanged unknowns d\ and d<± for unknowns a and b. 
What we gain from spin-flavor SU(6) is relations between the octet-octet coefficients q and 
the octet-decuplet coefficients dj. A determination of a subset of the coefficients q from 
low-energy octet-octet scattering would permit an estimate of the desired octet-decuplet 
coefficients dj. 

In summary, we have calculated the leading-order shift in the self-energy of the spin-| 
decuplet baryons in nuclear matter. Our work differs in two ways from earlier calcula- 
tions of the A isomultiplet self-energies in nuclear matter; we use chiral SU(3) symmetry 
to extend the calculation to include £*, H*, and fi~baryons and find new momentum- 
independent contributions from four-baryon operators in the effective Lagrangian. We 
have identified quantities independent of the coefficients dj, the helicity-splitting of the A 
and S* self-energy shifts and the resonance width of the A isomultiplet, which are pre- 



sented in Figures |3.2| and R73l In section |3.4j we discuss the origin and resolution of 



a logarithmic divergence of the A self-energy near nuclear saturation density, which is 
particularly relevant for future work when the coefficients d\ and d% are known. The two 
major short-comings of our results are that the coefficients dj of the four-baryon operators 
have not yet been determined and possible effects of the large nucleon-nucleon scattering 
lengths have not been included. 
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Chapter 4 

Heavy Kaon/Eta Effective Theory 

In this chapter an effective field theory for treating kaon and eta interactions in a non- 
relativistic framework is developed. The relevant effective Lagrangians are derived to 
an order sufficient for one-loop calculations of irrj, irK, and KK scattering processes. 
Coefficients in the effective Lagrangians at low orders are determined from a matching 
calculation with SU(3) and used to predict KK scattering phase shifts. 

In section |4 . 1| we present the motivation for developing the non-relativistic theory, 
describe the three key ideas which form the foundation, and outline the overall program 
to be followed. Sections |4.2| and [4.3| establish the elements and principles from which the 
effective Lagrangian is developed then detail the construction of the first several orders 
of the Lagrangian. The coefficients of the effective Lagrangian in the lowest orders of the 



chiral expansion are determined from a matching calculation in section [4.4|. A prediction 



of the s-wave KK scattering phase shift and brief concluding remarks on the outlook for 



the heavy kaon/eta effective theory are contained in section 4.5. 



4.1 Prospectus 

For the range of momenta Q < m n the rate of convergence of the chiral expansion in 
SU(3) chiral perturbation theory is limited by approximately the ratio of mass scales 
(m v /A x ) 2 ~ 0.3. For processes involving only pions with momenta Q < chiral pertura- 
bation theory based on a smaller SU(2)l x SU(2)r chiral symmetry provides an alternative 
description. One advantage of using SU(2) xPT for low-energy pion interactions is a better 
rate of convergence; generalizing the power counting arguments in section [Z72] the expan- 
sion parameter is approximately (m^/m^) 2 ~ 0.1. However, compared to SU(3) xPT, a 
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drawback of SU(2) xPT is the great reduction in the variety of processes to which the 
theory can be applied with only a slight reduction in the number of coefficients to be 
determined in the Lagrangian. Today, SU(2) %PT and SU(3) are considered com- 
plementary descriptions of pion interactions at low-energy and the relations between the 
coefficients in the two order-Q 4 effective Lagrangians are known |[20|| . 

Motivated by the comparison above, we seek an effective field theory tailored for low- 
energy interactions (Q < m^) of pions with kaons or eta mesons to exploit the relatively 
small ratio m^/rriK- One key ingredient of the effective field theory, as an expansion 
in 1/mjf, is a non-relativistic treatment of the kaon and eta degrees of freedom. Because 
the momentum scale dictates a relativistic treatment of the pion field, we cannot build the 
G x = SU(3)l x SU(3)r chiral symmetry of QCD into the effective Lagrangian explicitly; 
flavor rotations in SU(3)v would map relativistic pions into non-relativistic kaons and 
vice versa. However, we are able to include the smaller G' x = SU(2) L x SU(2) R chiral 
symmetry of QCD and a description of the pion degrees of freedom based on SU(2) x?T 
becomes the second key ingredient of the effective field theory. 

The non-relativistic field theory naturally divides into a number of n-body sectors 
distinguished by strangeness and the number of heavy fields n in the initial state. Sectors 
of the theory in which the number of heavy fields is not conserved present a problem 
for the momentum expansion; the annihilation of a heavy particle introduces pions with 
'hard' momenta (~m^ jr) ) in intermediate states and the power counting scheme breaks 
down. This problem afflicts sectors involving either KK pairs or the eta meson, which is 
subject to pion conversion as in T]K — > ttK. The exception to the general rule proscribing 
the eta meson is for a single eta meson interacting with any number of pions; in that case 
the eta meson is protected by G-parity ||55|| . (When isospin- violating effects are included, 
the decays rj — > 37r create new problems in the single-eta sector.) We focus on applying 
the effective field theory to the unaffected one- and two-body sectors, i.e., to iri], ttK, KK, 
and the sectors related by charge conjugation. 

When constructing the Lagrangian for the heavy kaon/eta effective theory, we must 
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introduce coefficients, or low-energy constants, which are not determined by the imposed 
symmetries. Even at low orders in the chiral expansion, several low-energy constants will 
appear in each of the independent sectors we consider. Due to the scarcity of experimen- 
tal data in some sectors, the standard method for extracting values for the low-energy 
constants is impractical. The alternative is to pursue a matching calculation onto a 
better-known theory at higher energy to determine the low-energy constants in terms of 
known parameters of the high-energy theory. In the regime of non-perturbative QCD, 
SU(3) xPT is the natural choice for the 'high-energy' theory. An additional benefit of 
matching onto SU(3) %PT is that the full G x chiral symmetry is implicitly restored to 
the effective Lagrangian through constraints on the low-energy constants. Because the 
matching calculation enhances the predictive power of the heavy kaon/eta effective the- 
ory, we consider it the third (and final) key ingredient of the theory. The essence of heavy 
kaon/eta effective theory is a reorganization of the chiral expansion in SU(3) x?T, keep- 
ing order- by-order only the terms relevant to the low-energy regime, which should result 
in improved convergence. 

The steps in the program we follow are to establish the effects of symmetry transfor- 
mations on the fields of the theory, detail the rules for constructing a minimal effective 
Lagranian, build the Lagrangian for the sectors of interest to order Q A , perform the match- 
ing calculation onto SU(3) x?T, extract values for the low-energy constants of the heavy 
kaon/eta effective theory, apply the results to scattering problems near threshold, and 
finally look for further problems to which the theory can be profitably applied. For the 
nK sector, this program has been started independently by Roessl |56| ; however, we find 
he has omitted terms from the effective Lagrangian at higher order (> Q 3 ). Several topics 
relevant to developing the effective field theory are addressed in the discussion: field redef- 
initions (or use of equations of motion); reparameterization invariance constraints on the 
effective lagrangian; reconciling relativistic and non-relativistic treatments in the match- 
ing calculations; and consequences of differences between the power counting schemes of 
SU(3) %PT and of the heavy kaon/eta effective theory. 
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To conclude a section titled "Prospectus" by considering future applications for the 
heavy kaon/eta effective theory seems appropriate, so we briefly skip ahead to the final 
step in the program. While the effective Lagrangians written out for the 7r.fr and KK sec- 
tors include isospin-violating terms for generality, the values of the associated low-energy 
constants have not been determined in the matching calculations. Extending the matching 
to determine the low-energy constants for isospin violation and electromagnetic interac- 
tions would be straight forward and suitable for studying processes such as irK — > 77r.fr 
and jK — > iruK near threshold. In another application, baryon isomultiplets from the 
spin-i octet and spin-| decuplet can be coupled to the heavy kaon effective theory with 
the coefficients in the Lagrangian matched to the G x -invariant theory in sections |2.3| , [2.4| . 
Unfortunately, in many interesting kaon-baryon systems, strangeness-exchange reactions 



do not conserve the number of heavy particles, e.g., NK — > S7r ^7|, resulting in a break- 
down of the power counting scheme; however, the NK and AK sectors are protected 
by strangeness conservation and are compatible with the heavy-kaon framework. On a 
more speculative note, phenomenological models inspired by and incorporating elements 
of the heavy kaon/eta effective theory may be developed for the problematic sectors 
like KK ||58|| . For example, similar work has been done on 7r + 7r~ atoms H59| , and by inte- 
grating out 7T° pairs and allowing for a non-hermitian effective Lagrangian |I0| . Finally, if 
experimental data for the threshold energy region become available in abundance, extract- 
ing the low-energy constants of the heavy kaon/eta effective theory from the data may 
provide insights into SU(3) x?T. Knowing values for the low-energy constants would per- 
mit verification of the constraints on the low-energy constants imposed by the matching 
calculation and would provide a means of testing the SU(3) chiral expansion. 
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4.2 Elements and Principles 

The Lagrangian of heavy kaon/eta effective theory is divided into sectors by the heavy 
fields present, 

c = c™ + c K + c* R + cr> + r KK + ■■■ (4.1) 

and the Lagrangian for each sector is subdivided as C x = £>f m a chiral expansion. All 
the purely pionic interactions are derived from C nn , where j — 2,4, ... , which is identically 
the Lagrangian of SU(2) %PT. The effective Lagrangians for interactions of a single heavy 
field with any number of pions are characterized by the chiral expansion j = 1,2,3, ... ; 
the expansion of £ nKK is similar except that the expansion starts with j = due to 
heavy-field contact terms. In this section we briefly establish our notation and identify 
the principles used to construct the effective Lagrangians. 

The building blocks of the heavy kaon/eta effective theory are introduced in direct 
analogy to SU(3) xPT. The three pseudo-Goldstone bosons of SU(2) %PT are intro- 
duced in the exponential representation; the most elementary constituents are u = e % ®l 2F 
and x = 2BM. q where 





^ V2n- -tt° J \ m d J 

and the parameters F and B are non-trivially related to the SU(3) %PT parameters F 



and Bq. (The relations between the parameters were obtained by Gasser and Leutwyler [2D 
we review them in section |4.4| .) In the purely pionic sector, the effective Lagrangian is 
written in terms of the fields U = u 2 and \ as in section |27|. Because the kaon fields are 
treated as matter fields which transform in the fundamental representation of SU(2)v, for 
coupling pions to kaons we adopt the notation used in section and reintroduce 



V, = l^d^u + ud^), (4.3) 
\ = i{u\u-vd u ^), (4.4) 
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X± = u] XU ] ± u~x}u. (4.5) 

Throughout this chapter this notation is reserved for the SU(2) versions of the fields. The 
heavy fields are the t] meson and the kaons, which form two SU(2)v doublets 



(4.6) 



As non-relativistic fields, the creation and annihilation operators are treated separately, 
thus we distinguish from K c . 

The symmetries which we explicitly build into the effective Lagrangian are G ' chiral 

A 

symmetry, parity (P), charge conjugation (C), and reparameterization invariance. The 
implementation of reparameterization invariance is discussed separately in the next sub- 
section. The transformations of the fields under G' , P, C, and hermitian conjugation (f) 





( K+\ 




( K° \ 






|, *.= ( 






\ K ° ) 




1 -K- / 



are collected for reference in Table ^J] along with the power of Q, i.e., the chiral index, 
associated with each field. Under each of the transformations listed, a covariant deriva- 
tive of a field DO transforms in the same way as O provided, in the case of parity, the 
Lorentz index of the derivative is also raised or lowered. The number of low-energy con- 
stants required is restricted sector-by-sector by imposing the constraints of hermiticity 
and invariance under G' x and P on the effective Lagrangian. Within sectors of non-zero 
strangeness charge conjugation does not restrict the form of the effective Lagrangian, but 
determines C ^ from £ K and C ™ from L^ KK . 

4.2.1 Reparameterization Invariance 

By treating the kaon and eta mesons as non-relativistic (heavy) fields we have performed 
an implicit rephasing, as in section |2.4j , compared to the relativistic counterparts K r 
and Tj r , 



K{x) ~ V2Me tMv - x K r (x), r]{x) ~ y/2M e m ' v - x K r (x). (4.7) 
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Table 4.1: Transformations of fields and under G' chiral symmetry, parity P, charge 
conjugation C, and hermitian conjugation f, and the chiral index associated with each 
field. 



Reparameterization invariance (RPI) is a consequence of arbitrariness in 

the choice of the velocity used in the rephasing. The choice of a different velocity in 
the relations between the relativistic and heavy fields is compensated by a shift in the 
residual momentum of the heavy field. Lagrangians written in terms of heavy fields 
rephased with different velocities must still give the same result for physical S'-matrix 
elements, which implies invariance of the Lagrangian under shifts of v^. A consequence of 
RPI symmetry in the effective Lagrangian is that coefficients of terms at different orders 



in the Q-expansion are related |62] and by building reparameterization invariance into an 
effective Lagrangian we reduce the number of low-energy constants to be determined. In 
this subsection we use the kaon as a representative example; all of the arguments apply 
equally well to the case of the eta meson. 

To construct a RPI Lagrangian it is necessary and sufficient, as shown by Luke and 
Manohar p2" |, that the rephasing velocity and derivatives of the kaon fields appear only 
in the following combinations: 



f -iMv» + DJK, 
[+iMv„ + D^K\ 



(4.8) 
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where we will refer to the operator V as the RPI-covariant derivative. V is a valid 
derivative in the sense that it satisfies the chain rule, 

DJ^OK) = (V^K j )OK + K*(DO)K + rfO(VK), (4.9) 

which permits the use of partial integration and total derivative arguments with respect 
to V . Only after building the effective Lagrangian in terms of the RPI-covariant derivative 
and expanding each term in powers of the mass M do we make any particular choice of 
frame, i.e., explicitly setting = (1,0,0,0). 

Superficially, developing the effective Lagrangian in terms of the RPI-covariant deriva- 
tive presents a problem for power counting. An infinite tower of possible terms, such 
as K^K, K'V V^ K , K'V V^V^V^K, . . . , will contribute to each order of the chiral expan- 
sion of the Lagrangian. On closer inspection, the fact that RPI guarantees a particular 
relationship between coefficients of terms in the effective Lagrangian assures that the 
effect of adding another term with more factors of V can always be compensated by a re- 
definition of the coefficients already present. As an explicit example we consider the effect 
of adding the term 5C = b KH^V V^V^K, where £ M represents an unspecified combination 
of light degrees of freedom, to the interaction Lagrangian £ int = a KH^V K. Expanding 
the combined result gives 

C int + 5JC = -iM{a-bM 2 )K\vi)K + (a - bM 2 ) KU^D^K (4.10) 
- 2bM 2 K?(v-e)(vDK) - ibM K^(y-£)D D^K + • • • , 

where we have used t> M t> M = 1. A subsequent change of a, a — > a + bM 2 , absorbs the con- 
tribution of 5C to the terms already present in £ int . The addition of another term, 
SC = cV x KH^V fj V u V u V X K , would be followed by a change of a and b to restore the form 
of the terms present before the addition of 5C . Generally, RPI makes the possibility of 
such redefinitions certain to all orders in the expansion of the effective Lagrangian and in 
every sector of the theory. 
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The procedure just proposed is an iteration of adding new terms to the effective 
Lagrangian then redefining the coefficients of the terms already present. An equivalent 
and more practical procedure is to subtract the contribution to terms already present in 
the effective Lagrangian from each new term when it is added, for instance 



RPI-covariant derivatives contracted with Lorentz indices on light degrees of freedom are 
never dropped. We implicitly make such subtractions from every term when we include 
the term in the effective Lagrangian. 

Having resolved the issue of infinitely many contributions to each order of the chi- 
ral expansion of the effective Lagrangian, we turn now to determining the proper power 
of Q to associate with the leading contribution of a term in RPI form. Occurrences 
of RPI-covariant derivatives are divided into two classes, those contracted with another 
RPI-covariant derivative and those contracted into the light degrees of freedom. In the ex- 
pansion of a contracted pair of RPI-covariant derivatives, the contribution from (±iMv) 2 
is cancelled by the subtracted terms; the leading contribution to the power counting 
comes from v-D so each contracted pair contributes Q l . The leading-order contribution 
from expanding a RPI-covariant derivative contracted with the light degrees of freedom 
is unaffected by the subtraction; the RPI-covariant derivative is counted as one power 
of or order Q°. To determine the effective Lagrangian to order Q D , in terms where the 
operators for the light degrees of freedom contribute d powers of Q and carry n Lorentz 
indices, we must consider some terms with as many as 



V^K ] V V V V V K 



(4.11) 




N = n + 2(D-d) 



(4.12) 



added RPI-covariant derivatives. Toward the upper limit of the necessary number of 
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RPI derivatives, some relief comes from the fact that total derivative arguments will 
allow replacing some terms with alternates where a derivative is moved onto the light 
degrees of freedom and the alternate terms may be neglected to the order we are working. 

In addition to the ambiguity in the choice of the rephasing velocity v^, we have free- 
dom to choose the mass M used in the rephasing relation, equation ( |4.7|) . This freedom 
was exploited in equation (|2.26|) to avoid awkward x^-dependent phases in the heavy 



baryon Lagrangian. The only constraint on the choice of M is that any residual mass SM 
is small to prevent failure of the power counting scheme. In the context of heavy 

kaon/eta theory, we perform the rephasing of the kaon and eta fields with the masses 
and M v respectively, where Mk and M v are defined as the meson masses in the SU(2) chi- 
ral limit m u ^ — > and m s finite. With this choice, the residual masses of the K and rj 
mesons at tree-level are generated through interactions with the field x+ an d are included 
perturbatively. Because the residual masses vanish in the SU(2) chiral limit, by definition, 
the effective Lagrangian does not include any explicit residual mass term —5MK^K; in 
this sense our choice for the rephasing mass is the 'natural' choice. 

4.2.2 Use of Field Redefinitions 

Redefinitions of the fields of the theory can be used in two ways to potentially improve 
or simplify the effective Lagrangian. The first class of field redefinitions eliminates terms 
from the effective Lagrangian which are proportional to a classical equation of motion; 
this class of field redefinitions, combined with total derivative arguments, is a powerful 
method for reducing the number of unknown coefficients required. The technique is well- 



documented in the literature |65|, |67|, |68| for simplifications of effective Lagrangians 
and we have nothing to add aside from identifying the general form of the terms we 
eliminate by this method. Terms with an explicit factor of D^A^, VV^K, V V^T], or the 
hermitian conjugates can be combined with contributions from other terms in the effective 
Lagrangian to form complete equation of motion terms and are subsequently eliminated. 
The second class of field redefinitions allows the elimination of all time-like derivatives 
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of the heavy fields from the interaction Lagrangian |JT], |35], |69|. The principle is closely 
related to the elimination of equation of motion terms, and if all equation of motion terms 
have already been eliminated by redefinitions of the first class then further redefinitions 
from the second class will permit a re-expression of the effective Lagrangian but not 
the further elimination of any low-energy constants present. Instead, the motive for 
considering field redefinitions in the second class is to simplify the application of the 
effective Lagrangian by replacing time-like derivatives of the heavy fields with spatial 
derivatives. Spatial derivatives permit a more transparent power counting for Q in loop 
diagrams than is possible when time-like derivatives occur in vertices. 

Considering the free field theory of a (complex) heavy scalar is sufficient to illustrate 
the general method and highlight conclusions relevant to how the matching calculation 
between relativistic and non-relativistic theories is performed. The free-field Lagrangian 
for the heavy scalar, after rephasing as in equation ( J4.7[ ), becomes 



iv-d + 



V\ + (iv-df 
2m 



(4.13) 



where the spatial derivative Vj_ is defined by d^d^ = (v-d) 2 — V^_. After substituting a 



field redefinition <f) Q = <fti + 5<ft, we find the well-known result that 
the time-like derivatives from terms at order 1/m and gives 



in i 



IV- 



cancels 



iv-d H — 

2m 



2,(ivdf (iv-d)V] 



16m 2 



Am 2 



+ C(l/m 3 



(4.14) 



(A slightly different approach is taken in references [^, [55].) Substituting a second 
field redefinition 0i = 02 + dependence on time-like derivatives is canceled in C^ 2 
through order 1/m 2 by 5cj)' = (4V^_ + 3(iv-d) 2 )<j)/32m 2 . Repeating this process, time-like 
derivatives of 4> can be eliminated to any arbitrary order in an expansion in terms of (d/m). 

Our original contribution to this discussion is the observation that the infinite series 
of field redefinitions can be formally summed and a relatively simple expression for the 
field redefinition <p = !F[(j)] can be obtained. The approach is based on substituting a 
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series representation for J-[(j)\, and deriving recursion relations from the constraint that 
all time derivatives cancel from Here we present a brief plausibility argument for 
the inverse solution (p = JF _1 [0 O ]; the actual recursion relations generated are presented 
in Appendix [FJ with a short description of their solution. The first step is to complete the 
square for the time-like derivative appearing in equation ( 4.13Q ; 



- ±6* 
_ 2 m ^0 



(iv-d + m) 2 - m 2 + V 2 ± 
iv-d + m — \/m 2 — 



(4.15) 



IV 



d + m + y/m 2 - V 2 



Second, observing (-y/m 2 — — m) is the relativistic kinetic energy operator K, we 
recognize the first expression in brackets is the kernel of the desired Lagrangian for <p, 
£- = (^(iv-d — K)(j). The second term in brackets is {iv-d + K + 2m) which we divide 
between O and 4>\ and finally use to identify </>; 



-iv-d + K + 2m 



2m 



iv-d — 



if -(9 + K + 2m 



2m 



(4.16) 



iv-d + K + 2m) /2m 



The explicit form for the redefinition resolves some differences between the perturba- 
tive field theories of the heavy fields 0o an d 4>- Developing the field propagator from the 
complete free-field Lagrangian in equation ( |4.13| ) results in two poles in the propagator, 
at iv -d = ±(m + K). Only one pole is relevant for the non-relativistic treatment and the 
second pole is eliminated differently in the 0o and <f) theories. The second pole is com- 
pletely absent from the Lagrangian £-; the disappearance of the second pole is traced to 
a singularity of the field redefinition !F[(f)] which causes <p to vanish at the location of the 
second pole. In the 0o theory, the pole is removed by treating the operator (iv-d) 2 /2m 
in as a perturbative correction to the O two-point function; the second pole does 
not appear at any finite order. The second difference between the two non-relativistic 
theories is the normalization of one-particle states relative to the normalization of the rel- 
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ativistic one-particle state |33|, [35], which becomes relevant when we pursue the matching 
calculation between the non-relativistic heavy kaon/eta theory and fully relativistic %PT. 
We determine the relative normalization by evaluating = JF _1 [0 O ] at the location of 



the relevant pole, (iv-d) = K = E — m, which gives ^2m(f) = y2E 0q. Which factor is 
included in matching calculations between relativistic and non-relativistic theories de- 
pends on whether relativistic corrections to the propagator are included as insertions 
of (iv-d) 2 /2m in the 0o theory or as insertions of V^_/8m 3 + V f ]_/16m 5 + • • • in the <fi 
theory. 

As stated above, the motivation for eliminating time-like derivatives of heavy fields 
is to simplify the power counting of loop diagrams involving both light- and heavy-field 
propagators. However, the field redefinitions required to eliminate time-like derivatives 
are not reparameterization invariant. We consider the question open as to whether RPI 
can be reformulated as a principle of the redefined (f) scalar field theory. Without explicit 
RPI symmetry of the effective Lagrangian, the relations between coefficients of different 
orders in the Q-expansion are not guaranteed. As a consequence, we forego using further 
non-RPI field redefinitions and treat the heavy kaon and eta fields in analogy with the 
0o scalar field theory described in this subsection. 



4.3 The Effective Lagrangians 

The approach we take to developing the effective Lagrangians is based on the work of 
Fearing and Scherer extending £ x pt to order Q 6 . We begin by addressing two points 
where we use a different approach, specifically when performing the SU(2)v contractions 
and when using total derivative arguments to eliminate terms. After discussing those 
issues we give a brief summary of the overall procedure then present the explicit con- 
struction of the effective Lagrangians C nK , C nri , and C 7vKK . 

Given a product of field operators, Fearing and Scherer form G x -invariant contractions 
by taking traces of all possible permutations of the matrix fields. (Roessl applies the 



47 



same method for C by introducing a matrix field KK> for the kaons pB[ .) Then trace 
relations, such as 



Tr[ABC] + Tt[ACB] 



(4.17) 



Tr[AB]Tr[C] + Ti[CA}Ti[B] + Tr[SC]Tr[A] - Tr[A]Tr[S]Tr[C] 



for 2 x 2 matrices, can be used to eliminate some of the resulting terms in favor of others 
in the set. For the relatively simple case of SU(2)v, we prefer to form a minimal set of 
contractions by appealing to the algebra for addition of angular momenta. We prefer the 
second method because 1) it permits a clean separation between isospin-conserving and 
isospin-violating terms, 2) the connection between generated terms and physical processes 
is more meaningful and direct, and 3) the method is more familiar. 

A general matrix field O transforms under SU(2)v as © 1, where the component 
is identified as Tr[0]. To separate out the 1 component of a matrix field, we define the 
notation 



Since Tr[Aj and Tr[Vj both vanish, these fields transform in just the 1 representation. 
This notation is particularly helpful when applied to x±'i we find that Tr[x + ] and \- are 
proportional to m, but Tr[x_] and x+ are proportional to (m u — rrid) and give isospin- 
violating terms. 

As a simple illustration of the method, the product of fields KK> A X- transforms 
as !®|®l<g)(0©l) which contains three singlet representations; the resulting set of 
contracted forms is Tr[x^A ]K*K, Ty[x-}K^A^K, and K'\x-,A ]K. We choose to take 
K^[x-, A^]K instead of the equally suitable form K^X-A^K because the commutator 
projects out just the 1 component of the product (%_ x A ). Also, we can distribute 
any covariant derivatives which are present after performing all SU(2)v contractions since 
they do not change the SU(2)y transformation of the field on which they operate. (How- 
ever, antisymmetric expressions such as [A , A ] must be kept until derivatives have been 



= 0-\Ti[0]. 



(4.18) 
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distributed across the fields.) 

The second point to discuss is how to select terms which are to be eliminated by adding 
total derivatives to the effective Lagrangian. Clearly, the set of selected terms should 
contain as few terms proportional to an equation of motion (EoM terms) as possible, 
since such terms will be eliminated anyway by subsequent field redefinitions. Fearing 
and Scherer p5[ take an approach driven by intuition, and apply rules such as all terms 
where more than half of the derivatives act on a single field may be eliminated, and use 
Lorentz index exchange arguments to preserve as many EoM terms as possible from the 
remaining set. Because our power counting scheme permits many more derivatives at any 
order than appear in C xPT , the intuitive approach is vulnerable to questions of whether 
every possible total derivative was used to eliminate a term and whether the set of terms 
retained contains the most EoM terms possible. 

We present an algorithm which reduces the use of total derivative arguments to linear 
algebra and assures that no EoM terms are eliminated at the expense of keeping another 
term (and low-energy constant) in the effective Lagrangian. We view the set of m terms 
prior to any eliminations as the basis of a m-dimensional vector space Q, and the set 
of coefficients of those terms in the effective Lagrangian represents an arbitrary vector 
in Q. Also, each total derivative, when expanded by the chain rule, is represented by a 
vector in Q; we define us to be the n-dimensional sub-space of Vt spanned by the complete 
set of total derivatives. The key ideas are that the vector representing the most general 
effective Lagrangian is arbitrary up to the addition of any vector from uj and any set 
of m — n terms which spans the complementary sub-space Q/u is sufficient for the most 
general effective Lagrangian. 

We construct a matrix P from the column vectors of a complete set of total deriva- 
tives, not necessarily linearly independent. A complete set of total derivatives for terms 
containing k derivatives can always be generated by explicitly taking the derivatives of 
each possible term with k — 1 derivatives. For any matrix O we define the operation NS[0] 
which returns a matrix whose columns span the (right) null space of O. Then, trivially, the 
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columns of the (m — n) x m matrix NS[P T ] span the desired sub-space Q/u. A particular 
choice of m — n terms from the original set also spans Vt/uj if the matrix R constructed 
from the column vectors representing the terms satisfies the relation det[P T • NS[P T ]] 7^ 0. 
Thus we can guarantee a complete set of total derivatives and test any particular choice 
of terms for suitability as a basis of Q/u. The remaining hurdle is to establish a means for 
identifying which choices of terms contain the maximum number of EoM terms; because 
the number of acceptable sets grows combinatorially with the number of derivatives and 
fields in each term, exhaustive testing of the possible sets is impractical. 

By definition the columns of NS[0] represent the linear combinations of the columns 
of O which vanish. The key observation is if we construct a matrix Q from the rows 
of NS[P T ] which correspond to the EoM terms, then the columns of NS[Q] represent all 
independent linear combinations of the columns of NS[P T ] with vanishing projection onto 
the EoM terms. The result is the columns of the matrix NS[P T ]-NS[Q] span the sub- 
space of tt/u which is orthogonal to all of the EoM terms, and the number of columns is 
exactly the minimum number of non-EoM terms which must be included in the effective 
Lagrangian. 

One considerable advantage of the algorithm just described is that the computations, 
from developing the complete sets of terms and total derivatives to selecting a single op- 
timal basis of terms to include in the effective Lagrangian, are easily programmed with 
software for symbolic mathematics. We also note that for very large sets of terms, the 
problem can be broken into two smaller problems by dividing the full set of terms consid- 
ered {tj}i<j< m into disjoint sets {tj + t'A and {i(tj — it)} and applying the algorithm to 
each smaller set. Finally, if the full set of terms forms a hierarchy of 'better' terms and 
'lesser' terms (beyond the distinction between EoM and non-EoM terms), the steps for 
removing a maximal set of EoM terms from the basis of Q/u can be repeated for each 
level in the hierarchy to assure the minimum number of terms must be selected from each 
of the successive 'lesser' sets of terms. 



In summary, we list the steps of the procedure, based on Fearing and Scherer [55 
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which we apply to construct the minimal effective Lagrangians C 7tK and £ nKK order by 
order in Q. 

1. Identify all simple (uncontracted) products of fields and derivatives consistent with 
parity and of the proper power of Q. 

2. For each product, form a minimal set of forms contracted over SU(2)v indices in 
analogy to addition of angular momenta. 

3. For each contracted form, distribute the derivatives over the fields and contract 
Lorentz indices all possible ways, keeping in mind that derivatives of x± are unnec- 
essary and D A v and all multiple derivatives are implicitly symmetrized. 

4. For sets of terms with many derivatives, the power counting may indicate that some 
contractions of Lorentz indices yield terms higher order in Q than the order to which 
we are working; drop all such terms. 

5. Replace each pair of terms t and with the hermitian combinations (t + tf) and 
i(*-ft). 

6. For each set of terms, apply the total derivative algorithm to select a basis set 
containing the most EoM terms possible. 

7. Eliminate EoM terms, proportional to either D^A^, V^V^K, or V^V^K^ through use 
of implicit field redefinitions. 

8. Expand RPI-covariant derivatives of the kaon fields as V = ±.iMkv u + D t then set 
^ = (1,0,0,0). 

In practice, after identifying a product of fields in step (1.), we work out all the terms 
generated by adding covariant derivatives to the product of fields up to the highest relevant 
order in Q before moving on to the next product of fields. 
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4.3.1 The tiK and tttj Sectors 

In the interest of transparency, so that the process and results can be independently 
verified and because some knowledge of intermediate results is necessary to extend the 
Lagrangian to higher order, we present details on some of the steps in the construction 
of the effective Lagrangian. After the discussion related to each product of fields, we list 
the RPI forms of the terms which are kept, together with the coefficient associated with 
each. The particular coefficient we assign to a given term depends on the leading power 
of Q contributed by the term after the RPI-covariant derivatives are expanded. 

The only product of fields permitted at order Q° is K^K. Because covariant derivatives 
can be freely moved from one field to the other, sophisticated total derivative arguments 
are unnecessary regardless of the number of derivatives distributed on the kaon fields; 
in each order of the Q-expansion, only a single term of this form needs to be retained. 
For no derivatives or two derivatives, the terms we keep are K^K and V^K^V^K. By 
definition, the coefficient of the mass term is the square of the mass used in the heavy 
field rephasing, i.e., the mass of the kaon in the chiral limit m u ^ d — > 0. For four or more 
derivatives, we can always move all derivatives except a single contracted pair off the kaon 
field leaving only the EoM term which is subsequently eliminated. The only combination 
of terms generated from the product K^K (at any order) is 

VK^K - M 2 K {K^K). 

No product of fields at order Q 1 is consistent with RPI and parity invariance. At or- 
der Q 2 two products of field operators must be taken into account, KK^x+ and KK^AA v . 
For the field product KK^x+ we use the contracted forms 

1i[ X+ ]K ] K (cf), K^x + K (cl K ). 

Because x+ is effectively 'transparent' to covariant derivatives (since D^x+ is equivalent to 
alternate terms we consider elsewhere), we can freely move covariant derivatives between 
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the kaon fields just as with the product K^K. Consequently, for any number of derivatives 
distributed on this form, we can always reduce the set of terms to a single EoM term 
proportional to V^V^K. 

For the field product KK^ A^A y we choose the contracted forms TrL4 AJ.fft'.K" and 
K* [A , A u ]K, and up to order Q 4 in each case we can distribute up to six derivatives over 
the fields. First we consider the form Tr[A ij^ K; adding two derivatives, we must 
keep a total of seven terms of which five may be EoM terms; adding four derivatives, we 
must keep a total of nine terms at order Q 4 of which eight may be EoM terms; adding six 
derivatives, we must keep a total of four terms at order Q 4 of which three may be EoM 
terms. The terms we select after all eliminations are 

Tt[A^]K^K 
Tr^ A"]VK^V V K {cl K ), 
Tr[D A V D^ A V \K^ K (ejj*), 
l^[A^D v A X \V x K^V i y v K + h.c. (e^), 
Tr[D fl A I/ D X A K ]V^V u K^V x V K K (e\ K ). 

The second form K^[A ,A\K vanishes unless derivatives are added; adding two deriva- 
tives, we must keep three terms of which two may be EoM terms; adding four derivatives, 
we must keep six terms of which five may be EoM terms; adding six derivatives, we only 
need to keep two EoM terms. The terms we select after all eliminations are 

VK\A\A V \VK {d?\ 
V^K^A", D V A X \V V V X K + h.c. {d\ K ). 

The only product of fields consistent with RPI and parity invariance at order Q 3 is 
KK^ A x— The product results in three contracted forms, Tr[x-A ]K^K, K^[x-,A]K, 
and Tr[x-]K^A K, which up to order Q 4 can be combined with either one or three 
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derivatives. Because all three of the contracted forms are a product of the same four 
distinguishable fields (no exchange symmetry), the total derivative arguments are exactly 
the same in each case; with one derivative we must keep two terms of which one may 
be an EoM term; with three derivatives we may keep only two EoM terms. The three 
non-EoM terms which we elect to keep are 

Tr[x_^]im^ + h.c. (df), 
iK^[x-,A^K + h.c. (df), 
Tr[ X -]K^ A^K + h.c. (d% K ). 

At order Q 4 four products of field operators are taken into account: KK^x+X+i 
KK^X-X--, KK " <A fl A u X+, and KK^ A^A v A x A k . The first two products will only ap- 
pear without any additional derivatives, and are both contracted the same way under 
SU(2)v- When coupling the representations of the two matrix fields, only the symmetric 
combinations [(0 ®1)®(0© l)] s = 0©0©l s ©2 will appear resulting in three SU(2) V 
contractions for each product. The six terms generated in this fashion are 

Tr[ X+ ] 2 ^ ( eir)j Tr[x- X -]K^K (e#), 
Tt[ X+ ]K^x + K (ejf), Ti[x^x-K «), 
Ti[ X+ X + ]K^K (e£), Tt[x^] 2 K^K «). 

The product of fields KK^A^A u x+ permits six distinct SU(2) V contractions and un- 
der the exchange A^ <-> A y three are symmetric and three are antisymmetric. For the six 
contracted forms, the power counting permits the addition of two derivatives. The three 
antisymmetric forms are Tt[ X+ A^A u }K^ K , Tr[ X+ }K^[A^ AJK, and K^[ X+ ,[ A ^ A M 
which all vanish if no derivatives are included. In each of the three cases, when two 
derivatives are included, the symmetry allows only a single term of order Q A which can 
be eliminated by total derivative arguments. The three symmetric contracted forms are 
TV[ X+ ]Tr[A A v }K^K, Tr[AA]K^ x+ K, and K\A X+ A u + A vX+ A)K. With the addi- 
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tion of two derivatives to each of the symmetric forms, we are required to keep only a 
single non-EoM term. The three symmetric forms result in a total of six terms in the 
effective Lagrangian at order Q 4 , 

Tt[ X+ ]Tt[A^}K^K (ejf), Trlx+M^ A V ]V^V V K «), 

Tr[A^]K^x + K (elf), 1t[A» A v }VK^x + V u K (etf), 

tfU^+^tf (e£), V^t(^x + ^ + ^x + ^)V^ (elf). 

The final product of fields, KK^A A V A X A K , transforms as §<g)§<g)l<g)l<g)l<g)l 
which contains nine singlet representations. Three of the nine SU(2)v contractions cor- 
respond to Tr[A^AjTr[A A AJi ; r 1 'i ; r under distinct permutations of Lorentz indices; the 
remaining six contracted forms are distinct permutations among the Lorentz indices of 
Tr[A AJK^[A X ,A]K. In each of the two cases, the power counting permits distribut- 
ing up to four derivatives on the fields, but the power counting also prohibits any EoM 
terms in the set of terms which can contribute at order Q 4 . First we consider the second, 
commutator-type form. The contributions of order Q A from this term vanish if either 
no derivatives or four derivatives are distributed on the fields; in the case of adding two 
derivatives, a single term of order Q A must be kept. The term we choose to keep is 

iTr^]^^, AV U V X K + h.c. {el K ). 

For the double-trace contracted form, there are two ways to contract the Lorentz indices 
without adding derivatives of the fields. In addition, with two derivatives on the fields 
we must keep two terms at order Q A ; with four derivatives we must keep a single term 
at order Q 4 . The net contribution to the effective Lagrangian is five terms; we select the 
following terms for this set: 

Tr [A^A 1 ] Tr \A v A v \ K^K (e} K ), 
Tt[AA]Tt[A^A"]K^K (eg*), 
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Tr[A^}Tr[A" A A ]V V K%K (eg*), 
r Fr[A A v ]Tr[A l 'A x ]V v KtV x K (eg*), 



We collect the terms of C nK listed above and present the complete results for the 
chiral expansion C nK = ^ . £* K to order Q 4 . We expand the RPI-covariant derivatives 
to separate the terms which contribute with different powers of Q. We also explicitly 
symmetrize the multiple derivatives and derivatives of A^ wherever necessary and absorb 
some constant factors into the coefficients where convenient. Finally, all of the isospin- 
violating terms are collected together at the end of each equation, offset with parentheses. 

C\ K = iK^D Q K (4.19) 
M K Cl K = \D^D^K + cl K Tx[A°A {) }K^K + cl K Ti[A^A tl }K^K (4.20) 

+ cl K Tr[ X+ ]K*K ( + dfrfx+k) 
M 2 K Cl K = tcl K TT[A°A^](K^D^K -h.c.)+tdl K K^A ,D A°}K (4.21) 
+ idl K (K*[A°, A»\DK + h.c.) + idl K Tr[x-A°}K^K 
+ dl K K^[ X -, A°]K ( + idfTr[x-}KU°K^ 
M z K Cl K = cl K Tr[A' i A v ]D K*D V K (4.22) 
- \d\ K \K\\A\D^\ + [A , -D Q A M ] + [A ,D»A°])D^ K - h.c.} 
+ dfDK\A\ A V \DK - \dl K ^\ X -A%K^DK - h.c.) 
+ {dt^Xx-.A^K - h.c.) +el K Ti[D A°D A°}K^K 
+ eZ K Tr[A' i D A ](KiDK + h.c.) + el K Tr[A A°]Tr[A A°]K^K 
+ el K Ti[D^A v (D^A u + D l 'A")]K^K + e^Tr^A^Tr^A^K 
+ eg* Tr [A°AJ Tr [A° A M ] K^K + e^Tr Tr [A v A v \ K^K 

+ el K Tr[A ll A v \Tr[A ,l A v \rfK + ie; K Tr[A° A^[A° , A»]K 
+ el«Tr[x + ]Tr[A A ]K^K + e^Trlx+MAA^K 
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+ elZK*A°x + A°K + egrfAJc+A»K + e^Tr[A°A°}K^ x+ K 
+ elfTr[A^]K^x+K + e^Tr[x + ]K^ x+ K + e^Tr[ X -}K^ X -K 
+ e^Ti[ X+X+ ]K^K + <Tr[ X _] 2 i^ 



This procedure does not need to be repeated for the nrj sector. The construction 
of £ 7r7? is identical to the construction of C kK with three exceptions: the eta transforms 
as a singlet under SU(2) V ; all isospin-violating terms must be omitted as discussed in 



section |4.1| ; and the effective Lagrangian within the irrj sector is invariant under charge 
conjugation. From the results for C nK we construct C nTJ directly by simply dropping 
all terms which violate isospin conservation or are not proportional to the general form 
(D m K^D n K), replacing M K by and each occurrence of (D m K^D n K) by (d m rfd n r)), 
and finally checking the remaining terms for charge conjugation invariance. 

= irfd r] (4.23) 
M V CZ> = ^d^d^rj + cY'Tr^A^T] (4.24) 
+ c^TrfAA"]^ + c^Trfx+l^ 

M*a? = icl v Tr[A A*](rfd^]-h.c.) + d nv Tr[ X -A }rfr] (4.25) 
Mffi = c^Tr^A^d^d^- |^Tr[x-A M ] (77^77 -h.c.) (4.26) 
+ el v Ti[D A D A°}^7] + e^TrfA" D Q A°}(rf 6^77 + h.c.) 
+ e^Tr^A^A" + D v A ti )]r^r] + e^"Tr[A c A°}Tt[A g <A°]rf rj 
+ e^>Tr[A A°]Tr[A A^ri + el v Tr[A A^]Tr[A°A^}rfr] 
+ e?Tr[A A^TrlAyA'tfTi + e\ T ' Tr [A A v ]Ti [A* 1 A u ]rfr] 
+ eSTr[ X+ ]Tr[^°A>^ + e^Tr [ X+ )Tr [A ^A^r] 
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4.3.2 The tiKK Sector 

In the tiKK sector, the full expression for the effective Lagrangian to order Q A is not 
necessary to treat 2 <-> 2 scattering processes until very high order. For instance terms 
involving four factors of only appear with at least two closed pion loops and will be 
suppressed by at least a factor of Q 6 . We develop the full effective Lagrangian to order Q 2 ; 
at orders Q 3 and Q 4 we consider only those terms which give rise to four-kaon contact 
terms. To distinguish contact terms from pion interactions, we note that an expansion 
of quantities in the number of pion fields gives ~ —d^n/2F, x~ — —^imBn/ F, and 
Tr[X-] ~ -U{mu - m d )B7i°/F. 

Also, we choose to couple the kaon annihilation operators together as KeK and 
KerjK, where e = ir 2 , before contracting with the other fields present for the advantage 
of greater symmetry under exchanges, antisymmetric and symmetric respectively. Cou- 
pling four kaon operators together in this way gives six representations under SU(2)v- 
The three representations involving an antisymmetric contraction of fields are simple 
products of the two forms above and transform as © 1 © 1; the remaining three 
combinations of fields are (KerjKy (KeTjK) (a 0), i€j kl (KeTjKy (Ker k K) (a 1), and 
((K eTj Ky (KeTkK) + (Ker^Ky (KeTj K)) (a 2), where use of e^, the three-dimensional 
Levi-Civita symbol, is always clear from context. For brevity of notation we define the 
combinations 

(K.K.lKsK,) = (K^K^iK^K,), (4.27) 
(^K^OIK.K,) = ie^K^K^TrirkOjiKsenK,) (4.28) 

where Kj represent general fields in the | representation of SU(2) V and O is any SU(2) V 
matrix operator. 

Without any factors of A^ or x± there are two contracted forms of four kaon field 
operators, (KeKy(KeK) and (K K\K K). To determine all contact operators up to 
order Q 4 requires distributing up to eight RPI-covariant derivatives on each form. (Many 
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with several derivatives contribute only at much higher order in tree diagrams; however, 
we include them for completeness regardless.) The algorithm for making optimal use of 
total derivative arguments serves us very well in this endeavor; we need to deal with as 
many as 221 candidate terms restricted by 179 linearly independent total derivatives. We 
consider the symmetric contracted form first and count the terms left after total derivative 
arguments; one term with no derivatives is necessary; adding two derivatives, we must 
keep three terms of which two are EoM terms; adding four derivatives, we must keep nine 
terms of which seven are EoM terms; adding six derivatives, 20 terms are required of 
which 18 are EoM terms; finally adding eight derivatives, 42 terms are required of which 
39 are EoM terms. After all redundant terms are eliminated, we select the following set 
of terms: 

(K K\K K) {a* KK ), (K V^K\K V^K) {b" KK ), 

iyKV il K\VKV v K) {q KK ), (yKVK\V ll KV v K) {q KK ), 

(VKV v V x K\V t KV''V x K) (dl KK ), (V,,KVV x K\V x KV' t V''K) [d\ KK ), 
(V^V U KV X V K K\V^V U KV X V K K) {e\ KK ), 
{V^V X KV V V X K\V K V^KV K V U K) {e\ KK ), 
(V^V U KV^V U K\ V X V K K V X V K K) (el KK ) . 

For the antisymmetric contracted form the term with no derivatives vanishes; adding two 
derivatives, we must keep a single non-EoM term; adding four derivatives, we must keep 
four terms of which three are EoM terms; adding six derivatives, twelve terms must be 
retained of which ten are EoM terms; adding eight derivatives, 28 terms must be retained 
of which 26 are EoM terms. After all redundant terms are eliminated, we select the 
following set of terms: 

(KeV^K) f (KeV^K) (c^ K ), 
{KtVK)\V u KeV^V u )+\i.c. (dl KK ), 
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(y u K6V^V x Ky(y X KeV^V v K) (dl KK ), 
(V v KeVjrK)t(V x KeV f *V x K) (eg™), 
(V^V u KeV x V K K) t ( V M V u KeV X V K K) (e\ KK ), 

{y u Kevy u Ky{y ) y R KeV ll v x v K K) . 



The final term in the above list does not contribute until order Q 5 , so we do not list an 
associated coefficient for that term. 

The next product of fields we consider is KKK^ K^x+ which results in five contracted 
forms and each form must be considered with up to four derivatives. Because x+ is trans- 
parent to covariant derivatives, three of the contracted forms yield terms analogous to a 
product of fields KKK^Kl The contracted forms Tr[x+](^ K\K K) and (K K\ X+ \K K) 
give four terms each after all redundant terms are eliminated in analogy to (K K\K K); 



The contracted form Ti\x+](KeK)\KeK) gives two terms after all redundant terms are 
eliminated in analogy to (KeK)\KeK)] 



Again, the final term in the list contributes only beyond order Q 4 . The final two contracted 
forms for the product KKK^K^x+ are the hermitian conjugate pair (KeKy (Ke X+ K) and 
(Kex+Ky (KeK) . Each of the forms vanishes without derivatives; adding two derivatives, 
in each case only one EoM term is required; adding four derivatives, four EoM terms are 
required. These final two contracted pairs do not result in any necessary terms in the 



Tt[ X+ ](K K\K K) (cf ™), 
Trl X+ ](KV^K\KV»K) (dl KK ), 
Tv[ X+ ]{VKV"K\V v KV v K) (eg™), 
Tr[ X+ ]{VKVK\V»KV v K) 



(KK\ X+ \KK) (cm, 
{KVK^KV^K) (dl KK ), 
{VKV»K\x + \VKV u K) (efD, 
{VKVK\ X+ \V»KV V K) 



Tr\x + }(KeV^(KeV»K) (eg™), 
Trix+^KeVKy^KeV^K). 
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effective Lagrangian. 

The only remaining terms we will consider at order Q A are derived from the field 
product KKK^ K^x+X+ without any derivatives. Considering the SU(2) V couplings gives 
eight distinct contractions, but four vanish due to antisymmetry. The remaining four 
terms contributing to the effective Lagrangian at order Q A are 

Tr[ X+ ] 2 (K K\K K) «*), Tr[ X+ ](K K\ X+ \K K) «*), 
Tr[ X+X+ ](K K\K K) «*), {Ke X+ K)\Ke X+ K) 

The final product of fields in the irKK sector is KKK^K^ A^A v and to derive the La- 
grangian to order Q 2 we only need to consider adding up to two derivatives. The SU(2) V 
contractions yield six independent forms, two of which are related by hermitian conju- 
gation. Without adding any derivatives to the product of fields, four of the contracted 
forms vanish due to an antisymmetry under an exchange of fields. Adding two derivatives 
to each of these four contracted forms gives terms which only appear at higher order 
than Q 2 ; we list the terms derived with the power of Q at which it would first contribute. 
From the hermitian-conjugate pair of contracted forms we get two EoM terms and 

(KeV^Ky(Ke[A»,A u ]V u K) (~Q 3 , twice), 
{KeV^iKe^D^A^K) (~Q 4 , twice); 

each form contributing both as (O + O^) and as i(0 — O^). The other two contracted 
forms with a pair of fields antisymmetric under exchange give two EoM terms and the 
higher-order terms 

{KVK\A»A V \KVK) (~Q 3 ), 
^\AA%KeVK)\KeVK) (~Q 4 ), 
Tx\A» A v \{KeVK)\KeV v K) (~Q 4 ). 
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The two terms without any derivatives which do contribute to C^ K are 

TilA^A^K K\K K) (c^), (KeA^K)^ (KeA^K) 

When two derivatives are distributed on either form, we must retain eleven terms of which 
five may be EoM terms. Of the twelve terms appearing in the effective Lagrangian, the 
only four contributing to C^ K are 

TrlA^A^KV^KlKV^K) (c\ KK ), 
Tr[A^A v }(K K\VK VK) + h.c. {cl KK ), 
{KeA^K^iyKeA v V v K) + h.c. (c} KK ), 
{KeA»VK)\KeA v V v K) + {KeA»V v K)\KeA v VK) {c\ KK ). 

From the product of fields KKK^K^A A v , the remaining eight terms derived are 

iTt\A^ A U ](K K\V K VK) + h.c. (~Q 3 ), 
i[KeA il K^{yKeA v VK) + h.c. (~Q 3 ), 
Tr[A^](K K\V U KV U K) (~Q 3 , twice), 
{KeA^iy^eA^VK) (~Q 3 , twice), 
^\DAD»A V \{KK\KK) (~Q 4 ), 
(KeD li A v K)^KeirA v K) (~Q 4 ). 

We collect the terms listed above and present the results for Cq KK through L\ KK . 
The contact terms of C\ KK are easily derived from the enumerated terms above; we do 
not present the fully-expanded form of C\ KK . 

M 2 K Cl KK = a nKK (K K\K K) (4.29) 
M\L\ KK = ib nKK {(K K\K D Q K) — h.c.} (4.30) 
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M A K Cl KK = V KK {K D^K\K D' I K) + cl KK (KeD fi K)\KeD >1 K) (4.31) 
+ {c\ KK + 2c\ KK ){K D K\K D K) 
- c\ KK {(K K\D Q K D Q K) + h.c.} 

+ (cl KK + cI kk )Tt[A q A q ](K K\K K) + cf K Ti[A^}(K K\K K) 
+ + c ; KK )(KeA K)\KeA K) + cl KK {KeAK)\KeA»K) 
+ cl* K Tr[ X+ }(K K\K K) ( + c\« K {K K\ X+ \K K)) 
M\Cl KK = icl KK {(K D Q K\D^K D^K) -h.c.} (4.32) 
+ 2ic\ KK {(K D^K\D K D^K) - h.c.} 

+ i(d\ KK + <£^ K ) | (K K\D Q K D D Q K) - (K D Q K\K D Q D Q K) 

-2(K D K\D K D K) - h.c. } 
+ idf K {(KeD Ky(KeD Q D K) - h.c.} 
+ id\ KK | {KeDK)\KeD Q D»K) + {KeDK)\KeD^D^K) 

+ 2{KtDK)\D Q KeD ii K) - h.c. } 
+ idl KK Tr\x + ]{(K K\K D Q K) - h.c.} 
( + idl KK {(K K\x + \K D Q K) - h.c.}) 
+ other terms involving A^ 

Having worked out the forms of the effective Lagrangians, we turn in the next section to a 
matching calculation to determine the coefficients which appear in terms of the parameters 
of the 'high-energy' theory, SU(3) %PT. 

4.4 Matching Calculations 

A topic which needs to be addressed before setting up the matching calculations be- 
tween the heavy kaon/eta theory and SU(3) chiral perturbation theory is the difference 
between the expansions of the two theories. The heavy kaon/eta effective theory is an 
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expansion in Q/M, where Q represents the generic scale of and external momenta \p |, 
and M represents the scale of heavy masses M K ^. In contrast, SU(3) is an ex- 

pansion in (Q/A) 2 where Q represents the scale of external momenta and the full set of 
pseudo-Goldstone boson masses. To reconcile the two effective field theories in the low- 
energy regime \p\ < m n , we must recognize three distinct mass scales Q, M, and A which 
represent respectively external momenta and m^, and M v , and the chiral symmetry 
breaking scale. 

Distinguishing the scales in is trivial. We compute scattering amplitudes to a 
given order in 1/A, then in each order separate the expression by powers of m s B ~ M 2 . 
In contrast, the role of A in the heavy kaon/eta effective theory is buried in the coefficients 
of the effective theory, just as the coefficients of are implicitly unknown functions 



of Aqcd and the heavy quark masses [13|. From within the heavy kaon/eta theory there 
is no way to determine the dependence of the coefficients on A via scaling arguments as in 
section pT2| because the dependence is obscured by the intermediate scale M. The matching 
calculation provides a direct way to establish the relationship between the coefficients and 
the chiral symmetry breaking scale. 

We match the theories by equating on-shell scattering amplitudes and the locations 
and residues of poles in the heavy field propagators of the two theories. (Off the mass- 
shell, the scattering amplitudes depend, in an unphysical way, on how the meson fields 
are defined and matching would require including all of the EoM terms eliminated by field 
redefinitions.) We choose to match amplitudes with a relativistic normalization; so a factor 
of \J 2Mk,t] is included for each external heavy-particle state in the heavy kaon/eta theory, 
consistent with the discussion in subsection |4.2.2| . For the amplitudes calculated in %PT 
we perform a non-relativistic expansion, making the replacement p° — > (Mk. v + k ). The 
expansions for an arbitrary amplitude A of mass-dimension d, in and in heavy 

kaon/eta theory respectively, are (symbolically) 

A xPT = M^ l ,^^ = M^^r i ^^ <"3) 

j,k j k 



64 



•Anr = Md Y, a ^V ( 4 - 34 ) 



where the (l/A)-expansion of xPT has been reorganized as an expansion in M/A within an 
expansion in Q/M. Equating the expansions in powers of Q determines the combinations 
of terms ajQi = J2 k aj^Q^ (M 2 / A 2 ) k for even j and (XjQi = for odd j. Solving for 
the coefficients in the heavy kaon/eta effective theory, we find that matching gives the 
coefficients as an expansion in powers of (M/A) 2 . This behavior is seen explicitly in the 
results for matching F and B to the parameters of SU(3) by comparing amplitudes 
for pion interactions. The matching relations were derived by Gasser and Leutwyler |20[ ; 



)SF^f + ^f + S(W}, (4.35) 

M 2 M 2 Ml ] 

B = M 1-^^1*^-16(^-2^)^ + C(M 4 /A 4 )j. (4.36) 

(x) 

The leading-order results from xPT determine coefficents of the effective Lagrangians Cj <3 
to leading order in (M/A) 2 (with the exception of coefficients for terms with time-like 
derivatives). Next-to- leading order results from xPT determine the same coefficients at 
the next order in (M/A) 2 and the leading order of a further set of coefficients from the 
heavy kaon/eta theory. 

We present the results for the leading order matching calculation then apply the results 
to make a prediction for KK scattering phase shifts in the next section. The parameters 
of C nK are, up to corrections of order (M/A) 2 , 

M 2 K = m s B , cl K = 0, cf' = -l/4, = -1/16, (4.37) 



and d\ K = d\ K = d\ K = 0. The parameters appearing in £ 7rr? are, up to corrections of 
order (M/A) 2 , 

M 2 = |m sJ B , cr = C2" = 0, C3" = -1/16. (4.38) 
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The coefficients d^ K and d KV are not determined until matching at next-to-leading order 
in xPT. For C wKK we only perform the matching for the coefficients of the four-kaon 
contact terms and find the leading order results are 

iykk _ - m sB KK _ —m s B K _ K _ m s B 

with corrections of order (M/A) 4 . 

In effect, the heavy kaon/eta effective theory becomes a dual expansion in Q/M 
and (M/A) 2 through the matching calculation. Applying the theory to physical processes 
requires establishing a relative weight to the two expansion parameters. The relative 
weight dictates to what order of the expansion in powers of (M/A) 2 one must work for 
a consistent result to a particular order in Q/M. Equivalently we introduce a common 
expansion parameter £ and assign Q/M and M/A each a characteristic power of £. The 



ratio of masses m^/Mx ~ 0.3 and the typical size of SU(3) symmetry breaking effects pO 
suggest assigning the ratios of mass scales Q/M and (M/A) 2 equal powers of the param- 
eter £, so we take £ ~ Q/M ~ (M/A) 2 ~ 0.3. 

4.5 Application: KK Scattering 

From the effective Lagrangians C kK and C wKK we calculate the KK scattering amplitudes 
in the isospin 1=0,1 channels to order Q 2 . Using the results of the matching calculations 
of section |4.4|, we make a prediction for the KK s-wave scattering phase shift in the 1=1 



channel. We plot a comparison of the leading order results of the heavy kaon/eta theory 
and SU(3) %PT and comment on the usefulness of this approach. 

Initially we calculate the KK scattering amplitude as an expansion in powers of Q/M, 
treating all coefficients as intrinsically C?(£°). The leading-order contribution to the scat- 
tering amplitude is order Q° and arises from the tree diagram with an cr Kif -vertex. Be- 
cause time derivatives on external legs contribute Q 2 , the effective Lagrangian results in 
no contribution from tree diagrams at order Q. At order Q 2 we include tree diagrams with 





Figure 4.1: Feynman diagrams contributing to the KK scattering amplitude in the heavy 
kaon/eta effective theory; (a) order Q° tree diagram, (b) order Q kaon-bubble diagram, 
(c) order Q 2 tree, double-bubble, and pion exchange diagrams. 

a vertex from C\ K or L\ KK and which give one of the coefficients \f KK , c% KK , or c\q k . 
Because the two-body sector admits nearly-infrared-divergent behavior, the appropriate 
power counting for loop diagrams is given in equation ( |2.31| ). The power counting shows 
that kaon-bubble diagrams are suppressed by a single power of Q, so to order Q 2 we 
take double- and single-bubble diagrams with only the a* KK -vertex included. Finally, 
one-pion loop diagrams are suppressed by only the powers of Q generated at the vertices, 
which means double pion exchange is included with vertices from the order-Q effective 
Lagrangian C\ . The set of relevant diagrams is illustrated in Figure |i~T . 

Integrations over loop energies dq are performed by contour integration and the result- 
ing momentum integrals dq are evaluated in dimensional regularization. The scattering 
amplitudes are calculated in the center-of-mass frame and the order-Q 2 equation of mo- 
tion, 2MxEk = k 2 — 16c^ K rhB, gives the kaon energy of the external states. The results 
for Aff R are 

iA$ R {k) = 32icf«^-, (4.40) 

M K 



iA%{k) = 32ia- { 1 + ia-A - } (4.41) 

„„ Eh r^rhB rriiMl- m 2 

+ mV KK — + 2hQic\n % w K In —1. 

M K M 2 K 16tt 2 F 4 /i 2 
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From the matching calculation we determined c^ KK ~ 0; consequently, the (1=0) p-wave 
scattering amplitude and phase shift vanish to the order we are working. Including the 
^-scaling of the coefficients, the leading contribution to A$ R comes from the order-Q° tree 
diagram and counts as order £. The next terms in the phenomenological expansion in £, 
of equation (}4.41|) , are the order-Q kaon-bubble diagram and the order-Q 2 tree diagram. 
To work consistently at this order would require extending the matching of the coefficient 
in the leading-order result, a 7vKK , by two orders in (M/A) 2 . 

For comparison, the result for the scattering amplitude in x?T is 



l * PT ~ 3F 2 



o 



8M% + 12k 2 - 2{m s + m)B . (4.42) 



The leading-order results for the s-wave phase shift in both the heavy kaon/eta theory 



and SU(3) are presented together in Figure fO|. The s-wave KK scattering length 



determined from the phase shift plotted in the figure is 

a« = 0.45 x lCr 13 cm, (4.43) 



which is consistent with a repulsive K-K interaction potential. In both cases the leading 
corrections are suppressed by ~(M/A) 2 ~ 0.3. In light of the fact that the heavy kaon/eta 
theory is determined directly from the scattering amplitudes of x?T, the agreement of 
the two expansions to within the expected 30% corrections is not surprising. 

The calculation of the s-wave phase shift in both xPT and the heavy kaon/eta theory 
has illustrated a fundamental point relevant to the utility of the matching calculation. 
The motivation for developing the heavy kaon/eta theory was to achieve a better expan- 
sion in the threshold regime by virtue of the better ratio of scales (Q/M) 2 ~ .1 versus 
(M/A) 2 ~ .3. In effect, the price for the improved convergence was an increase in the num- 
ber of low-energy constants, i.e., coefficients in the Lagrangian, to be determined. The 
initial proposal was made to determine the unknowns of the low-energy theory by perform- 
ing a matching calculation onto a theory with fewer unknown parameters, SU(3) %PT. 
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k McV- 1 




(b) 



ft McV - 

Figure 4.2: (a) Leading-order calculations of the KA' s-wave (1=1) scattering phase 
shift So(k) in degrees as a function of center-of-mass momentum k; the dashed line is 
heavy kaon/eta theory, the solid line is SU(3) x?T. (b) Relative error between the results 

K R -sf T \/\sn 

However, by matching onto we guarantee that the heavy kaon/eta theory is limited 
to converging no better than the theory to which it is matched. Indeed, we find that 
the matching calculation explicitly reintroduces the mass-scale ratio (M/A) 2 that we at- 
tempted to circumvent. The matching calculation appears to be both the savior and the 
Achilles' heel of the heavy kaon/eta effective theory. 

For a wide variety of processes, including KK scattering, the heavy kaon/eta effective 
theory does not provide a computational advantage over xPT, and the additional work 
required to carry out the matching calculation cannot be justified. However, the outlook 
for the heavy kaon/eta theory or for the use of matching calculations is not entirely bleak. 



The approach is useful is two senses that were suggested in section |X The first is that 
independent of the matching calculation, when data are available for a direct empirical 
determination of the low-energy constants, then taking advantage of the better rate of 
convergence in the non-relativistic theory will be possible. For a short-term perspective, 
the available data on tcK scattering [170, [7JJ may be useful in this way. 

The second way in which the non-relativistic effective theory may be useful is calculat- 
ing quantities for which the relativistic formalism of x?T is awkward. An example where 
a non-relativistic treatment has proven beneficial is calculating the properties of bound- 
state systems || |32|, |7^, |73]|. In applications of this sort, reviving the matching calculation 
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approach has merit because the motive for using the non-relativistic effective theory is 
not solely to improve the convergence of the expansion. Since bound state formation 
always involves summing an infinite number of Feynman diagrams, the power counting 
considered here cannot be directly applicable for these problems. Two cases should be 
considered regarding bound states in the heavy kaon/eta effective theory. First, the effec- 
tive field theory cannot be applied to deeply bound systems because they lie outside the 
range of the momentum expansion. Second, for shallow bound states, the power counting 
scheme needs to reflect the fine-tuning which is implicit in the associated large scattering 
length. In systems with a shallow bound state, the matching calculation will not provide 
an estimate of coefficients of the terms responsible for binding the system, but may per- 
mit the determination of coefficients which contribute perturbatively to properties of the 
bound state. The interesting possibilities for strong-interaction bound states, KK and 
KN, both suffer kaon annihilation and are not suitable candidates. Another possibility 
is to calculate strong-interaction perturbations to a Coulomb bound state with a kaon, 
for instance tt~K + or K~H + , or kaon electromagnetic form-factors via e~K + . 

In conclusion, we have presented the foundations and laid out the principles for an 
effective field theory to describe the interactions of pions with non-relativistic kaons and 
eta mesons. The effective Lagrangians for few-body sectors were constructed explicitly 
for the first several orders in a chiral expansion. Much work remains to be done both 
in applying the theory to systems of interest and in generalizing the theory to include 
electromagnetic interactions, couplings to baryons, and the study of isospin violation. 
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Appendix A 

Symbolic Expansion of jC xPT 



This appendix presents Mathematicd^ routines 74] for the symbolic expansion of the 



SU(3) chiral lagrangian in terms of the embedded chiral SU(2) fields in the sincere 
hope that they will prove useful to others. The source code is divided into two files, 



lagrangian. math on page [72] which specifies the sequence of steps in expanding and 
saving the desired parts of the lagrangian and definitions .math on pages |73|-f?4| which 
defines the functions called to perform the expansion and simplify the results. Currently, 
the chiral Lagrangian C 2 in equation (|2.8|) is expanded out to terms including six boson 
fields; £ 4 in equation (|2.10p is expanded to four boson fields. The source code can be 



applied to higher-order parts of the chiral Lagrangian or to terms involving more bosons 
by adding transformation rules to the trace routine in definitions .math at line 42. 

The routines are of very limited use without a brief description of the symbols appear- 
ing in the input and output files. In addition to the obvious parameters F, B, and L1-L8 
defined as in section pT2l the quark masses are included as 



mhat — > m = \{m u + m d ), mdif f — > {m u — m d ), ms — > m s 
The correspondences for the boson fields are 

pion — > 7i t , kaon — > K, kbar — > K, eta — > rj. 



Each occurrence of a boson field is multiplied by a tag, boson, which simplifies expand- 
ing the lagrangian and separating the lagrangian into parts by the number of bosons 
1 Mathematica is a registered trademark of Wolfram Research, Inc. 
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appearing. Also, two matrices are given symbolic names in the input files; one refers to 
the 2x2 identity matrix and tau3 refers to the Pauli matrix r 3 = a z , which appears in 
isospin violating terms. Finally, the following notation appears in the output: 

• derivatives, d[A,x] — > d x A, 

d[pion,mu] — > (d fl 7r a )r a , d[kaon,nu] — > d u K 

• isovector products, v[A,B] — > ^Tr[AB], 

v [pion, pion] — > a 7i b Ti[r a r b ] = n ■ ir, v[d[pion,mu] ,tau3] — > d^n° 

• isospinor contractions, m[A,B,C] — > A a B a pCp, 
m[kbar,d[kaon,mu]] — > Kd^K, m[kbar ,pion,kaon] — > 7r ■ (KtK) 

• commutators, c[A,B] — > [A, S], 

v [c [pion, d [pion, mu]] ,tau3] — > 7r o (^7r 6 )(2ie ab3 ) = 2i (7r x (9^7?) • r 3 , 



m 



[kbar,c [d[pion,mu] ,d[pion,nu]] ,kaon] — > 2i (9 7r x c^7r) • (KtK) 



The output is written to the file lagrangian. save and breaks the lagrangian into 
pieces of manageable size. The different parts of the lagrangian are named LXyZ, where 
X is the chiral order as in C 2 versus £ 4 , Z is the number of boson fields appearing in the 
terms, and y is a letter which indicates the number of factors of M. q occurring in place 
of derivatives. A few select lines of lagrangian. save are given below to illustrate. 

L2a2 = d[eta, mu] "2/2 + m[d[kbar, mu] , d[kaon, mu] ] + 
v[d[pion, mu] , d[pion, mu]]/2 

L2b2 = eta~2*(-(B*mhat)/3 - (2*B*ms)/3) + (-(B*mhat) - B*ms)*m[kbar , kaon] - 
B*mhat*v [pion, pion] 

The translation of the lines above is 

Aa2 = \ ■ d»n + d^KcfK + \ d^ V , 

£-2b2 = —Bom n ■ n — Bo(m s + m) KK — |£> (2m s + m) rj 2 . 
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(******************** lagrangian . math ************************) 
(* symbolically expands the SU(3) ChPT lagrangian in terms of the *) 
(* embedded SU(2) fields, out to 6 bosons for L2 and 4 bosons for L4 *) 
« "definitions .math" ; timetag ["Session started :"] ; 

mdiff=0; (* drops all isospin violating terms, remove line to keep them *) 
X=Xd=2*B*{{mhat*one+mdif f *tau3/2 , 0} , {0 , ms}} ; 

P=boson*Sqrt [2] *{{pion/Sqrt [2] +eta*one/Sqrt [6] , kaon} , {kbar , -2*eta/Sqrt [6] }} ; 
U=matrixseries [Exp [(I/F) *#] &,P, {boson, 6}] ; Um=derivative [U,mu] ; 
Ud=matrixseries [Exp [(-I/F) *#] &,P , {boson, 6}] ; Udm=derivative [Ud,mu] ; 
timetag ["Defined the basics:"] ; 

(* Calculate L2a *) 

UmUdm=product[dot,Um,Udm,{boson,6>] ; L2a=trace [UmUdm] *(F/2) ~2; 

{L2a2 , L2a4 , L2a6}=Map [reduce [Coefficient [L2a, boson ,#]]&, {2 , 4 , 6}] ; 

L2aN=Expand [reduce [L2a/ . boson->l] - (L2a2+L2a4+L2a6) ] ; 

Save [ " lagrangian . save " , L2a2 , L2a4 , L2a6 , L2aN] ; 

{Um , Udm , UmUdm}=Map [truncate [# , {boson , 4}] & , {Urn , Udm , UmUdm}] ; 

Clear [L2a,L2a2,L2a4,L2a6,L2aN] ; timetag ["Calculated L2a:"]; 

(* Calculate L2b *) 

XUd=dot[X,Ud] ; UXd=dot [U,Xd] ; L2b=trace [XUd+UXd] * (F/2)~2 ; 

{L2b0 , L2b2 , L2b4 , L2b6}=Map [reduce [Coefficient [L2b , boson , #] ] & , {0 , 2 , 4 , 6}] ; 

L2bN=Expand [reduce [L2b/ . boson->l] - (L2b0+L2b2+L2b4+L2b6) ] ; 

Save [ " lagrangian . save " , L2b0 , L2b2 , L2b4 , L2b6 , L2bN] ; 

{XUd , UXd}=Map [truncate [# , {boson , 4}] & , {XUd , UXd}] ; 

Clear [U,Ud,L2b,L2b0,L2b2,L2b4,L2b6,L2bN] ; timetag ["Calculated L2b:"] ; 

(* Calculate L4a *) TUmUdm=trace [UmUdm] ; 

TUmUdn=trace [product [dot.Um, (Udm/.mu->nu) , {boson, 4}]] ; 

L4a=(Ll*product [Times ,TUmUdm, (TUmUdm/ .mu->nu) , {boson, 4>] 
+L2*product [Times , TUmUdn , TUmUdn , {boson , 4}] 
+L3*trace [product [dot , UmUdm , (UmUdm/ . mu->nu) , {boson , 4}] ] ) ; 

L4a4=reduce [Coefficient [L4a, boson, 4] ] ; L4aN=Expand [reduce [L4a/ .boson->l] -L4a4] ; 

Save ["lagrangian. save" ,L4a4,L4aN] ; 

Clear [Um, Udm, TUmUdn, L4a,L4a4,L4aN] ; timetag ["Calculated L4a: "] ; 
(* Calculate L4b *) TpXU=trace [XUd+UXd] ; 
L4b= (L4*product [Times , TUmUdm , TpXU , {boson , 4}] 

+L5*trace [product [dot, UmUdm, (XUd+UXd) , {boson, 4}]] ) ; 
{L4b2 , L4b4}=Map [reduce [Coefficient [L4b , boson ,#]]&, {2 , 4}] ; 
L4bN=Expand [reduce [L4b/ . boson->l] - (L4b2+L4b4) ] ; 
Save [ " lagrangian . save " , L4b2 , L4b4 , L4bN] ; 

Clear [UmUdm, TUmUdm, L4b , L4b2 , L4b4 , L4bN] ; timetag ["Calculated L4b: "] ; 
(* Calculate L4c *) TmXU=trace [XUd-UXd] ; 

TpXUXU=trace [product [dot , XUd , XUd , {boson , 4}] +product [dot , UXd , UXd , {boson ,4}] ] ; 
L4c= (L6*product [Times , TpXU , TpXU , {boson , 4>] 

+L7*product [Times, TmXU.TmXU, {boson, 4}] +L8*TpXUXU ) ; 
{L4c0,L4c2,L4c4} = Map [reduce [Coefficient [L4c , boson, #] ] &,{0 ,2 ,4}] ; 
L4cN = Expand[reduce[L4c/.boson->l]-(L4c0+L4c2+L4c4)] ; 

Save ["lagrangian. save " ,L4c0,L4c2,L4c4,L4cN] ; timetag ["Calculated L4c:"] ; 
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definitions .math ***********************#) 
(* defines functions called by "lagrangian .math" *) 

(* time-stamp printing function *) 

timetag [msg_] := (Print [] ; timenew=SessionTime [] ; 

If[untimed, (timeorigin=timeold=timenew; untimed=False) ; ]; 

Print [ StringTake [msgOspaces , 20] ,PaddedForm[(timenew-timeold) ,{9,3>] , 

PaddedForm[(timenew-timeorigin) ,{9,3}] ,PaddedForm[TimeUsed[] ,{9,3}] ] ; 

timeold=timenew ) ; 
untimed=True ; spaces=" "; 

(* catagories of symbols which are used in simplifying rules *) 
f ields={pion,kaon,kbar ,eta} ; ident={{one , 0} , {0 , 1}} ; 
structs={pion,kaon,kbar,one,tau3}; bras={kbar}; kets={kaon}; 

(* functions called directly from "lagrangian. math" *) 

truncate [exp_ , xpn_] : =Expand [Normal [Series [exp , {First [xpn] , , Last [xpn] }] ] ] ; 
derivative [exp_ , var_] : =(D [exp, var ,NonConstants->f ields] 

//. Dot[___,0,___]->0 //. Literal [D] [f_,v_, ]->d[f,v] ); 

dot [x_ , y_] : = (distribute [{ 

{Dot [x[[l,l]] ,y[[l,l]]]+Dot [x[[l,2]] ,y[[2,l]]] , 

Dot [x[[l,l]] ,y[[l,2]]]+x[[l,2]]*y[[2,2]]}, 
{Dot [x[[2,l]] ,y[[l,l]]]+x[[2,2]]*y[[2,l]] , 

Dot [x[ [2,1]] ,y[[l,2]]]+x[[2,2]]*y[[2,2]]} }] ) ; 
product [mult_ ,x_ ,y_ ,xpn_] := (Module [{i , j , coef f ,result=0} , ( 
coeff=Coefficient[#l, First [xpn] ,#2]&; 
For [i=0 , i<=Last [xpn] , i++ , For [ j =0 , j <=i , j ++ , 

result +=mult [coef f [x, j] , coef f [y , (i-j )] ] *First [xpn] ~i ; ] ; ] ; 
Return [Expand [result] ] ) ; ] ) ; 
matrixseries [func_,matr_,xpn_] := (Module [{x,i=0,m={truncate [matr,xpn] }}, ( 
While [Union [Flatten [Last [m] ] ] ! ={0} , 

(AppendTo [m, product [dot , Last [m] , Fir st [m] , xpn] ] ; i++) ; ]; 
Return [ReleaseHold [Expand [ truncate [func [x] ,{x,0,i}] 

+Limit [f unc [x] ,x->0] * (Hold [ident] -1) ] /. x~n_ . : >Hold[m[ [n] ] ] ]]);]); 
trace [exp_] : = (Expand [ExpandAll [tr [exp [[1,1]]] +tr [exp [[2,2]]]] 
// . {tr [x_? (FreeQ [# , any [structs] ] &) ] ->x , 

tr [x: (_Plus I _Times)] : >Map [tr ,x] , tr [x_~n_Integer] ->tr [x] ~n } 
/ . tr [x_Dot] : >isoscalars [Apply [List ,x] ] 

//.{tr[one]->2, tr [tau3 I pion I d [pion, _] ] ->0 , tr [x_ ,y_] ->2*v [x,y] , 

tr [x_ , a ,x_] ->v [x,x] *tr [a] , tr [a ,x_ ,x_ ,b ] ->v [x,x] *tr [a,b] , 

m[a ,x_,x_,b ] ->v [x,x] *m [a,b] } 

/ . {tr [x_ , y_ , z_] ->v [c [x , y] , z] , tr [x_ , y_ , x_ , y_ , z_] ->2*v [x , y] *v [c [x , y] , z] , 
tr [x_ ,y_ ,z_ ,w_] ->2*v [x,y] *v [z ,w] -2*v [x,z] *v [y ,w] +2*v [x,w] *v [y ,z] , 
m [a_ , x_ , y_ , b_] ->v [x , y] *m [a , b] +m [a , c [x , y] , b] II , 
m[a_,x_,y_,z_,b_]->v[c [x,y] ,z] *m[a,b] /2 

+v [x , y] *m [a , z , b] -v [x , z] *m [a , y , b] +v [y , z] *m [a , x , b] , 
m [a_ , x_ , y_ , x_ , y_ , b_] ->v [x , y] *m [a , c [x , y] , b] 

+2*v [x,y] *v [x,y] *m[a,b] -v[x,x] *v[y ,y] *m[a,b] } 
/. {v [tau3 ,tau3] ->1 , m [_ , c [x_ ,x_] , _] ->0, 

v [c [x_ , y_] , z_] : >0/ ; ( (x===y ) I I (x===z) I I (y===z) ) }] ) ; 
reduce [exp_] : = (Module [{terms , forms , f reeq=FreeQ [# , any [f ields] ] &} , ( 
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terms=Map [commonf orm, Flatten [{ 

Replace [ExpandAll [exp] ,x_Plus : >Apply [List ,x] ] }] ] ; 
forms=Union [DeleteCases [terms , _?freeq, {1}] //. _?f reeq*x_->x] ; 
Return [Collect [Apply [Plus , terms] .forms] ] ) ; ] ) ; 

(* functions used internally to simplify and reduce expressions *) 
any[list_] : =Apply [Alternatives , list] ; 

indexed [form_] :=(MatchQ[form, ( any[structs] I d [any [structs] , ] 

I Dot [_? (FreeQ [#, any [bras] ]&) ,__] I Dot [__,_? (FreeQ [#, any [kets] ] &)] )]) ; 

distribute [exp_] := (ExpandAll [exp] //. Dot [ ,0, ] ->0 

//.{Dot [w ,x_Plus,y ] : distribute [Dot [w,x,y] ,Plus] , 

Dot [w , a_* (x_?indexed) , y ] ->a*Dot [w , x , y] , 

Dot [w ,one,y ]->Dot[w,y] }) ; 

isoscalars [list_] := (Module [{olist , break=Map [First .Position [list , any [bras] ] ] } , ( 
If [break==-Q, Return [Apply [tr, list]] ;] ; 
olist=RotateLef t [list .First [break] -1] ; 
break=Append [break-First [break] +1 .Length [olist] +1] ; 
Return [Product [ Apply [m, Take [olist , {break [ [i] ] .break [ [i+1] ] -1}] ] , 

{i.l, Length [break] -1} ]]);]); 
commonf orm [exp_] := (Module [{vars , swaps .array} , ( 

array={l , exp} // . List [a_ , b_? (FreeQ [# , any [fields] ] &) *x_] ->{a*b , x} 

//.{{a ,x_*y_,z }->{a,x,y,z}, {a ,x_~n_Integer ,y }->{a,x~ (n-1) ,x,y}}; 

array=Apply [Union [Flatten [Outer [Times , ##] ] ] &, Map [Union, Map [List , array] 
/. {{m[a_,c[x_,y_] ,b_] }->{m [a, c [x,y] ,b] , -m[a, c [y ,x] ,b] } , 

{v [c [x_ , y_] , z_] }->{v [c [x , y] , z] , -v [c [y , x] , z] , 

v[c[y,z] ,x] ,-v[c[z,y] ,x] ,v[c [z,x] ,y] ,-v[c [x,z] ,y] }, 

{v [x_ , y_] }->{v [x , y] , v [y , x] } } 
/ . {m [a_ , c [tau3 ,x_] ,b_] ->-m [a, c [x,tau3] ,b] , v [tau3,x_] ->v [x, tau3] , 

(v[c[tau3,x_] ,y_] I v [c [y_ ,tau3] ,x_] )->v[c[x,y] ,tau3] }] ] ; 
vars=Union [Cases [exp, (d [_ ,x_] ->x) , Infinity] ] ; 
swaps=Map [Thread [Rule [vars ,#]]&, Permutations [vars] ,{1}] ; 
Return [First [Union [Flatten [Map[(#/. swaps) &, array]]]]]) ; ]) ; 
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Appendix B 

Recursion Relations for (j) — > T\(\) 



We start with the free-field Lagrangian for a (complex) heavy scalar field, 



iv-d + 



V 2 ± + (iv-df 
2m 



where is the scalar-field velocity and d = (v-d) 2 — V^_. We seek a redefinition of the 
scalar field = which eliminates secondary time-like derivatives of the heavy scalar 
in favor of higher powers of the space-like derivative Vj_, 



TV 



7 = 1 x ' 



6. 



The general field redefinition we consider here is 



j=0 k=0 



iv-d 



ml \m 



72 \ k 



and constraints on the form of give recursion relations for the coefficients f k . Using 
integration by parts, we find the recursion relations are 



( /0n2 _ i fl__lfO fO _ _ 1 fl 

\Jo) ~ l i JO ~ 4^0) Ji ~ 2^0? 



f fj+1 _ lfOfjil / fP f j-p i 9 fP+1 fj-p\ 
JOJO ~ 4JO JO ~r 4 / , Uoip 'Jo JO )i 



p=0 
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-tiff 
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JoJk+i 



9k+l 



3-1 k 

— I P f _i_ I \ \ ( f 'P fj-p-l I o fP+1 fj-p-1 I fp+1 fj~P\ 

2J0Jk+l ' 4 / j / j Wq-'fc-q+l ' q ■> k—q+1 •> q J k-q) 
p=0 q=0 

3-1 k 

+ I (fk+lfo P + 2fk+lfo P ) + | ^ {fqfk-q fq+lfk-q) 

p=0 q=0 
k-l 

\{fofk + fofk) + I {fqfk-q + /g /fc-g + ^fq+lfk-q) ) 

1 i'O /•O 

2 / <? ■/ fc— qi 
q=0 



and determine the coefficients in the sequence /q , f$ j /i > /o > /i > if > /o > • • • 

Solving for .F[0] is simplified by recognizing that the coefficients gj must reproduce 
the kinetic energy and well-known relativistic corrections; C-r = ffiiiv-d — K)4> where 



K 



m* 
V 2 



vi 



m 



2m 8m 3 16m 5 

We postulate that the field redefinition is a function of one of the combinations (iv-d ± Vj_) 
or (iv-d ± K) only, then verify that only the form (iv-d + K) satisfies the recursion rela- 
tions. We construct the solution through trial and error, guided by comparing expansions 
of the trial function with the recursion relations above, and find 



1 + 



iv-d + K 
2m 



Inspired by hindsight we present a formal, and much shorter, derivation of this result in 
Chapter f| 
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